arXiv: 1503.03317v2 [astro-ph.CO] 16 Sep 2015 


Probing Dark Energy models with neutrons 


Guillaume Pignol 

LPSC, Universite Grenoble Alpes, 

CNRS/IN2P3, Grenoble, 

France* 

(Dated: September 17, 2015) 

There is a deep connection between cosmology - the science of the infinitely large -and 
particle physics - the science of the infinitely small. This connection is particularly 
manifest in neutron particle physics. Basic properties of the neutron - its Electric 
Dipole Moment and its lifetime - are intertwined with baryogenesis and nucleosynthesis 
in the early Universe. I will cover this topic in the first part, that will also serve as an 
introduction (or rather a quick recap) of neutron physics and Big Bang cosmology. Then, 
the rest of the manuscript will be devoted to a new idea: using neutrons to probe models 
of Dark Energy. In the second part, I will present the chameleon theory: a light scalar 
field accounting for the late accelerated expansion of the Universe, which interacts with 
matter in such a way that it does not mediate a fifth force between macroscopic bodies. 
However, neutrons can alleviate the chameleon mechanism and reveal the presence of the 
scalar field with properly designed experiments. In the third part, I will describe a recent 
experiment performed with a neutron interferometer at the Institut Laue Langevin that 
sets already interesting constraints on the chameleon theory. Last, the chameleon field 
can be probed by measuring the quantum states of neutrons bouncing over a mirror. In 
the fourth part I will present the status and prospects of the GRANIT experiment at 
the ILL. 
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I. THE INTERPLAY BETWEEN NEUTRON PARTICLE 
PHYSICS AND COSMOLOGY 

There is a profound connection between cosmology - 
the science of the infinitely large - and particle physics 
- the science of the infinitely small. The primordial cre¬ 
ation of the light nuclei is an example: a detailed knowl¬ 
edge of nuclear interactions is needed to predict the out¬ 
come of the Big Bang nucleosynthesis. Conversely, cos¬ 
mological observations have a great impact on our under¬ 
standing of the laws of Nature at the most fundamental 
level. The existence of Dark Matter, Dark Energy and 
the matter-antimatter asymmetry provide three glimpses 
of what is hiding beyond the Standard Model of parti¬ 
cle physics. It is likely that future progress in particle 
physics will proceed from the solution of one of these 
cosmological puzzles, or vice versa. 

The deep connection between particle physics and 
cosmology is particularly manifest in neutron particle 
physics, when we use low energy neutrons to investigate 
the fundamental interactions and symmetries. The work 
presented in the next chapters builds on this tradition. 
We will contemplate the possibility to probe models of 
Dark Energy with neutron experiments. Before devel¬ 
oping on this original idea, we review in the present 
chapter the traditional interplay between neutron par¬ 
ticle physics and cosmology. A comprehensive treatment 
of this topic can be found in the excellent review of Dub- 
bers and Schmidt (2011). See also Dubbers (2014) for a 
shorter and more recent overview. 

There are two fundamental properties of the neutron 
connected to the physics of the early Universe. First, the 
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value of the neutron lifetime is a basic input for the 
theory of Big Bang nucleosynthesis. Second, the search 
for the neutron electric dipole moment is key to elu¬ 
cidate the origin of the baryon asymmetry of the Uni¬ 
verse. We will explain these two subjects, after a brief 
review of neutron particle physics and cosmology. 

A. Free up neutrons 

Neutrons make up about half the mass of ordinary ob¬ 
jects, like bicycles, cheese or elephants. Even an empty 
bottle, that is, full of air, contains more than 10 20 neu¬ 
trons per cm 3 . These neutrons are bound into nuclei 
since billions of years. On the contrary free neutrons do 
not exist naturally because they decay with a radioactive 
period of ten minutes. It takes a significant effort to ex¬ 
tract neutrons from nuclei; it is much more difficult than 
extracting electrons from atoms. Indeed the binding en¬ 
ergy of neutrons in nuclei is of the order of 1 MeV whereas 
the chemical binding energy is only about a few electron- 
Volts. Thus, no chemical reaction would set neutrons 
free: a nuclear reaction is needed. 

Neutron sources. About twenty large nuclear in¬ 
stallations dedicated to neutron production are presently 
being operated over the world. These installations fall 
into two categories: fission nuclear reactors and spalla¬ 
tion sources. 

The reactor of the Institut Laue Langevin (ILL) in 
Grenoble belongs to the first category. Neutrons are 
produced by fission in a compact cylindrical core (con¬ 
taining 9 kg of weapons grade uranium, 93 % 235 U) and 
moderated in heavy water. Each fission ejects on aver¬ 
age 2.5 neutrons with an energy of 2 MeV. The reac¬ 
tor design optimizes the thermal neutron flux, reaching 
10 15 cm' 2 s _1 with a reactor thermal power of 58 MW. In 
a typical nuclear reactor used to produce electricity with 
a thermal power of 3000 MW, the neutron flux is about 
10 14 cm' 2 s _1 . Among the 246 operating research reac¬ 
tors listed by the International Atomic Energy Agency, 
46 of them are high flux reactors (above 10 14 neutrons 
cm -2 s _1 ), 11 of them feature a cold neutron source and 
operate instruments available to external users (see Table 

I)- 

The ultracold neutron facility at the Paul Scherrer In¬ 
stitute (PSI) in Switzerland uses a neutron source based 
on the spallation process. A proton beam, with an en¬ 
ergy of 600 MeV and an intensity of 2.5 mA is shot at a 
lead target. When a proton collides with a heavy nucleus 
(lead in this case), it ejects about 20 neutrons with an 
energy of about 20 MeV along with other nuclear frag¬ 
ments. The neutrons are then moderated in a heavy 
water tank surrounding the spallation target. Spallation 
sources are now competing with reactors as neutron fac¬ 
tories, Table II lists the spallation sources equipped with 
neutron instruments available to external users. 


TABLE I High flux reactors producing beams of cold neu¬ 
trons. The number of instruments available for outside users 
is indicated. 


Reactor 

City 

Started 

Th. Power 

Instruments 

ILL 

Grenoble (France) 

1971 

58 MW 

45 

HFIR 

Oak Ridge (USA) 

1965 

85 MW 

12 

CARR 

Bejing (China) 

2010 

60 MW 

12 

FRM II 

Munich (Germany) 

2004 

20 MW 

19 

HANARO 

Daejon (Korea) 

1995 

30 MW 

8 

WWR-M 

Gatchina (Russia) 

1959 

18 MW 

19 

NIST 

Gaithersburg (USA) 

1967 

20 MW 

28 

JRR-3M 

Tokai (Japan) 

1990 

20 MW 

26 

BRR 

Budapest (Hungary) 

1959 

10 MW 

11 

OPAL 

Sydney (Australia) 

2006 

20 MW 

7 

BER II 

Berlin (Germany) 

1973 

10 MW 

22 


TABLE II Neutron spallation sources serving instruments 
available to users. 


Spall. Source 

City 

Beam Power 

Instruments 

SINQ @PSI 

Villigen (Switzerland) 

1.5 MW 

20 

SNS @ORNL 

Oak Ridge (USA) 

1.4 MW 

25 

JSNS @KEK 

Tsukuba (Japan) 

0.3 MW 

17 

ISIS @RAL 

Oxford (UK) 

0.2 MW 

36 

Lujan ©LANSCE 

Los Alamos (USA) 

0.1 MW 

13 


What is the purpose of investing so much efforts and 
money for building sources of free neutrons? What is 
special about free neutrons that cannot be done with the 
neutrons inside nuclei? Beyond doubt the main desir¬ 
able feature of the neutron is its electrical neutrality 4 . 
Even atoms, which have no net electric charge, interact 
mainly with electromagnetic interactions due to their rel¬ 
atively large polarisability. Atomic processes are essen¬ 
tially unaffected by nuclear or gravitational interactions 
and all phenomena we experience in everyday life are ulti¬ 
mately due to basic electromagnetic interactions (except 
for the rather trivial effect of weight). Using neutrons one 
can probe matter using non electromagnetic interactions. 
Studying the diffraction of neutrons at a sample mate¬ 
rial, solid-state physicists are able to obtain information 
about the structure of the sample which is complemen¬ 
tary to the information they get using electromagnetic 
probes such as X-ray diffraction. 

Particle physicists also find neutron’s sensitivity to 
non-electromagnetic forces quite useful to do fundamen¬ 
tal physics. Studying free neutron beta decay we mea¬ 
sure the fundamental parameters of the weak interac¬ 
tion. Measuring neutron bounces we explore gravity with 
quantum objects. The search for the neutron electric 
dipole moment constitutes a crucial test of the time rever¬ 
sal invariance in elementary processes. One can also use 


1 The particle data group (Olive et al., 2014) evaluation of the 
neutron charge is (—2±8) X 10 -2J electron charges, in agreement 
with electrical neutrality. 
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neutrons to search for new, yet undiscovered interactions. 
All these experiments would be impossible with charged 
particles for which the Coulomb force overwhelms any 
other force by many orders of magnitude. We shall come 
back to these subjects in the rest of this manuscript. 

At the ILL, 38 instruments are dedicated to the study 
of material structure, condensed matter and magnetism 
and 7 instruments are dedicated to quantum, nuclear and 
particle physics. For a broad exploration of the physics 
of slow neutrons, see Byrne (1994). 

Wave particle duality. Experiments very often 
benefit from using slow neutrons. Those neutrons result 
from the thermalisation of the primary fast neutrons in 
a thermal or cold moderator. In a moderator, a neutron 
will dissipate its energy by collision until it reaches the 
kinetic energy (kT = 25 meV at room temperature) of 
the molecules. The De-Broglie wavelength A of a neutron 
with kinetic energy E is 


A = 


2ttH 

\j2mE 


= 0.18 nm x 


25 meV 


E 


( 1 ) 


where m = 939.6 MeV/c 2 is the neutron mass. It is 
remarkable that the wavelength of thermal neutrons cor¬ 
responds to the typical distance between atoms in solid 
matter. Thus, whereas fast neutrons behave essentially 
as particles, cold neutrons could behave like waves. 

E. Fermi was the first to realize that slow neutrons 
could undergo optical phenomena such as reflection and 
refraction at surfaces. It can be shown (see for exam¬ 
ple Golub et al. (1991)) that matter acts as a uniform 
medium with a potential energy for slow neutrons called 
the Fermi potential. The Fermi potential is given by the 
following expression: 


V F = 47t — bn, (2) 

2m 

where b is the bound coherent neutron scattering length 
of the nuclei constituting the material and n is its num¬ 
ber density. If the material is heterogeneous one must 
sum the potentials of all nuclear species composing the 
material. For most materials the Fermi potential is pos¬ 
itive, i.e. repulsive 2 ; it is of the order of 10“ 7 eV, much 
smaller than the kinetic energy of thermal neutrons. For 
example, natural nickel is a material with a relatively 
high Fermi potential which amounts to 250 neV. 

Neutrons approaching a surface at grazing incidence 
could be reflected by the Fermi potential of the material. 
Total reflection occurs when the kinetic energy associated 


2 This is surprising because the strong nuclear interaction respon¬ 

sible for the Fermi potential is attractive: it holds neutrons inside 

nuclei. For an explanation of this apparent paradox, see for ex¬ 

ample Pignol (2009). 


with the velocity normal to the surface is smaller than 
the potential barrier (2). In practice this phenomenon 
is at play in neutron guides. It is possible to transport 
neutrons from the core of the reactor where they are pro¬ 
duced to an experimental hall situated at a distance of 
up to 100 m, using evacuated rectangular tubes, with 
a cross section of typically 100 cm 2 , made up of plates 
coated with nickel or with a multilayer of nickel and tita¬ 
nium. Thermal neutrons with an energy of E = 25 meV 
are reflected off a nickel surface if the angle of incidence 9 
satisfies the condition sin 2 9 <V F /E, that is, 9 < 0.2 deg. 
For colder neutrons, the efficiency of the guiding is better 
and a beam of cold neutrons transported in a guide can 
have a larger angular divergency. 

Ultracold neutrons. Zel’dovich (1959) real¬ 

ized that neutrons with total kinetic energy lower than 
250 neV should be reflected at any angle of incidence 
and therefore could be stored in material bottles. These 
storable neutrons were called ultracold neutrons (UCNs) 
by the pioneer workers in the field. It is not easy to get 
ultracold neutrons, because their proportion in the en¬ 
ergy spectrum of thermal neutrons (when neutrons are 
thermalised in a moderator at room temperature) is only 
10 -11 . The first storage of ultracold neutrons was re¬ 
ported by Groshev et al. (1971). The group of Russian 
physicists stored on average 2 neutrons at a time in a 
copper cylindrical bottle (diameter 14 cm and length 
174 cm). The capacity of the bottle to store neutrons 
was quantified by a storage time of about 30 s (compare 
to the neutron beta decay lifetime of 880 s). 

The kinetic energy of ultracold neutrons corresponds 
to a temperature as low as a few mK. However, it is 
important to understand that ultracold neutrons can be 
stored in bottles at room temperature. The neutrons do 
not thermalise with the walls of the bottle, like photons 
of visible light do not thermalise when reflecting off a 
mirror. In fact the specular (i.e. mirror) reflection of 
ultracold neutrons and visible photons share many simi¬ 
larities, because the wavelength of UCNs is of the order 
of 100 nm, very close to the wavelength of visible light. In 
both cases the wavelength is much larger than the lattice 
spacing of atoms in the matter of the walls. It means that 
the particle interacts with a large number of atoms in the 
wall, it is almost blind to the thermal motion of individ¬ 
ual atoms. Using appropriate wall materials, UCNs can 
undergo as many as 10 4 * specular reflections before being 
inelastically scattered or captured by a single nucleus. 

The ability to store neutrons for a long period was 
immediately recognized as a great opportunity for fun¬ 
damental physics, for two main reasons. First, ultracold 
neutrons provide a direct method to measure the neu¬ 
tron lifetime by simply storing neutrons for a certain time 
and counting the survivors. Second, longer observation 
times enhance the sensitivity of detection schemes like 
the Rabi or Ramsey resonance techniques, following the 
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uncertainty relation SE x T > h, where SE is the pre¬ 
cision of the measured energy and T is the observation 
time. Aim long apparatus observes a neutron spin in 
a cold beam for typically 10 ms, whereas in a storage 
experiment observation times longer than 100 s could 
be achieved. The measurement of the neutron electric 
dipole moment (nEDM) could thus benefit from an im¬ 
provement of the sensitivity by four orders of magnitude. 
Since 1970 the techniques to produce and handle ultra¬ 
cold neutrons were greatly improved to pursue these two 
goals: a precise measurement of the neutron lifetime and 
the search for the nEDM. We will come back to these two 
topics after taking a cosmological detour. 


B. The standard picture of the early Universe 


The modern picture of the Big Bang was initiated by 
G. Lemaitre, supported by the measurement of the ex¬ 
pansion of the Universe by Hubble (1929). At that time, 
the modern theory of gravity, Einstein’s general relativ¬ 
ity, was just invented. Within this framework it was 
possible to consider the Universe as a physical system 
interacting with its content. Assuming the cosmologi¬ 
cal principle - the Universe is homogeneous and isotropic 
on large scales - the cosmic evolution of a spatially flat 
Universe is described by the Robert son-Walker metric 


ds 2 = dt 2 — a(t) ( dx 2 + dy 2 + dz 2 ) (3) 


where a(t) is the dimensionless scale factor defined in 
such a way that at the present time to the scale fac¬ 
tor is unity a (to) = 1. In our expanding Universe, 
the Hubble parameter H(t) = a(t)/a(t) is positive. In 
the cosmic past, all distances were smaller by the fac¬ 
tor a(t). For example, light emitted at time t with a 
wavelength A is observed now with a larger wavelength 
Ao = X/a(t) = {z + 1)A. The redshift z = — 1, being 

actually observable, is often preferred to t as a variable 
to indicate the ticking of the cosmic clock. 

Assuming that the Universe is filled with an homo¬ 
geneous fluid with an energy density p and a pressure 
p, Einstein’s equations of general relativity (with a van¬ 
ishing cosmological constant) reduce to the Friedmann 
equations 


H 2 


8t tG 


3 


P 


dp 

dt 


+ 3H{p + p) = 0. 


( 4 ) 

( 5 ) 


In the early Universe, the expansion was dominated 
by radiation, i.e. photons and relativistic particles. In 
this case the pressure is related to the energy density as 
p = p/3. Then from the Friedmann equations we deduce 
p oc a _4 (f), a(t) oc \Jt and H(t) = X-. Going backwards 
in time (higher z) the energy density of the radiation 
p oc (z + l) 4 was higher. Assuming that the radiation was 


in thermal equilibrium, one can associate a temperature 
T to the radiation according to 

P = ^ N ( T ) T ^ (6) 

where N{T) counts the effective number of relativistic de¬ 
grees of freedom (+1 per boson polarization state, +7/8 
per fermion state). Combining with the Friedmann equa¬ 
tions we can relate the age of the Universe t to the tem¬ 
perature T 


t = 


45 


16t t 3 GN(T) T 2 


2.4 s 1 MeV 2 

y/W) T2 


( 7 ) 


The Universe was hotter in the past. For example, at the 
beginning of the Big Bang nucleosynthesis the tempera¬ 
ture was T ss 1 MeV, the effective numbers of degrees 
of freedom was N(T) = 29/4 (only photons and three 
families of neutrinos), thus the age of the Universe was 
t ss 0.5 s. One can also relate the photon temperature to 
the redshift, 


T = {z + 1)T 0 


N(To) 

N(T) 


1/4 


( 8 ) 


where To = 2.7255(6) K is the present temperature of 
the Cosmic Microwave Background (CMB). Note that 
Eq. (7) is valid only in the radiation dominated epoch 
whereas Eq. (8) holds for the whole thermal history of 
the Universe. 

Very different energy scales of the primordial plasma 
were successively at play in the thermal history of the 
Universe. This leads to a fascinating and fruitful con¬ 
nection between cosmology and particle physics that re¬ 
late phenomena at the highest energies with phenomena 
at the earliest epoch. Before exploring more specifically 
the relevance of low energy neutron physics, which is a 
subfield of particle physics, we shall shortly review the 
present standard scenario of the hot Big Bang (see Ta¬ 
ble III) . For a more complete review on the status of the 
standard cosmological model, see Bartelmann (2010). 

It is believed that the Universe started with a phase of 
accelerated expansion driven by a scalar field called the 
inflaton (Guth, 1981). What happened before is even 
more speculative because it is related with physics at the 
Planck scale and involves the great problem of unifying 
gravity with quantum mechanics. The inflaton field is 
supposed to be coupled to Standard Model particles and 
the energy density of the scalar field would be converted 
into particles by the end of the inflation. This reheat¬ 
ing process would start the thermal history of the Uni¬ 
verse with an initial temperature of presumably about 
10 15 GeV. 

Then the Universe cooled down while expanding and 
a succession of phase transitions occurred. At tem¬ 
peratures corresponding to the electroweak scale (T ss 
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TABLE III The history of the Universe: standard scenario. Values in bold are taken from Olive et al. (2014). 




Transition 

Temperature Time 

Redshift 2 



_ Inflation ends (reheating) 

~ 10 15 GeV 

~ 10' 37 s 

~ 10 28 


$ 

Electroweak phase transition 

170 GeV 

10- 11 s 

10 15 



QCD phase transition 

170 MeV 

30 ps 

10 12 

- 

1 

Big Bang nucleosynthesis 

80 keV 

150 s 

3 x 10 8 



Matter radiation equality 

1 eV 

30000 Yr 

3360(70) 



CMB decoupling (recombination) 0.3 eV 

372000 Yr 

1090.2(7) 



First stars (reionization) 

30 K 

460 MYr 

11 



v\ Today 

2.7255(6) K 

13.81(5) GYr 0 


200 GeV), the Higgs boson field expectation value con¬ 
densed from the symmetric phase (</>) = 0 down to the 
non-symmetric phase (0) = 246 GeV. We will come back 
to the electroweak phase transition in section I.D. Next, 
at a temperature below T = 100 MeV, quarks and gluons 
condensed into nucleons during the QCD phase transi¬ 
tion. The Universe was then populated by photons, neu¬ 
trinos, relativistic electrons and positrons, and traces of 
protons and neutrons. The nucleons are leftovers from a 
tiny asymmetry between quarks and antiquarks, the so- 
called baryon asymmetry of the Universe (BAU) which 
was generated earlier. Further condensation of nucle¬ 
ons into nuclei became possible at the temperature of 
T k, 80 keV during the Big Bang nucleosynthesis (BBN). 
We will explain the connection between the BBN and the 
measurement of the neutron lifetime in the next section. 

The radiation dominated period of the Universe ended 
at a temperature of T = 1 eV, when the energy den¬ 
sity of dark matter became larger than that of radia¬ 
tion. Shortly after, at T = 0.3 eV, electrons and protons 
combined into hydrogen atoms and the Universe became 
transparent to photons. When we observe the CMB to¬ 
day we see in fact photons emitted during the recombi¬ 
nation at the redshift of 2 ~ 1000. In the matter domi¬ 
nated Universe, the primordial density fluctuations grew 
until star formation began at the very recent redshift of 
z = 11. In the late luminous Universe the expansion 
is again accelerating, perhaps driven by a new substance 
called Dark Energy. We will develop in the next chapters 
some possible ways to address the Dark Energy problem 
with neutron experiments. 


C. Big Bang nucleosynthesis and the neutron lifetime 

The Big Bang nucleosynthesis. Gamow (1948) 
formulated the hypothesis that nuclear fusion in the early 
Universe is responsible for the formation of nuclei. He 
argued that the primordial nucleosynthesis started at a 
temperature of Then = 100 keV, when the radiative dis¬ 
sociation of deuterium stopped. The age of the Universe 
was then a few minutes, which he considered to be the 


timescale associated with nucleosynthesis. Conveniently, 
this is also the order of magnitude of the neutron life¬ 
time, which was poorly known in 1948. Then, he re¬ 
quired the mean time for the reaction p + n —> d + "/ to 
be about 2 minutes. Indeed, if the mean time was much 
longer then no complex nuclei would have been formed 
during the BBN. On the other hand if the mean time was 
much shorter then there would be no deuterium left (sub¬ 
sequent fusion reactions happen to be faster than deu¬ 
terium formation). Quantitatively this educated guess, 
sometimes referred to as the Gamow criterion, can be 
expressed as 

rib av x 2 m i n ~ 1, (9) 

where a « 10” 29 cm 2 is the fusion cross section, v ss 
5 x 10 8 cm/s is the thermal velocity of the protons and 
neutrons at the temperature Teen and is the number 
density of protons and neutrons. With the above crite¬ 
rion Gamow got rib ~ 10 18 cm~ 3 corresponding to a mass 
density of p b (t BBN ) ~ 1CT 6 g/cm 3 . 

Following this line of argument one can pre¬ 
dict the present temperature of the CMB, T 0 = 
Teen \/Pb(t 0 )/pb(t BBN ) , by using the scaling laws T(t) a 
a(f) _1 and pb{t) oc a(t)~ 3 . Using the estimate pb{t 0 ) ss 
10~ 3 ° g/cm 3 obtained from the Hubble rate, Alpher and 
Herman (1949) predicted To ~ 5 K. The Big Bang model 
was unambiguously confirmed in 1965 when Penzias and 
Wilson discovered the CMB. 

Building on the seminal work of Gamow, a more so¬ 
phisticated theory of the BBN was worked out. A key 
quantity is the ratio of baryon to photon number densi¬ 
ties r) = rib/n-y. For a given value of p the theory can 
predict the relative abundance of the four light elements 
produced during the BBN, namely D, 3 He, 4 He and 7 Li. 
From the observed relative abundances of each of these 
elements one can extract p. The most sensitive probe of p 
comes from the measurement of the deuterium to proton 
ratio in extragalactic clouds, D/H= (2.53 ±0.04) x 10 —5 , 
from which the following result is obtained (Cooke et al ., 
2014) 


??BBN = (6.0 ± 0.1) x 10 10 . 


( 10 ) 
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The study of the temperature anisotropies in the CMB 
provides an independent determination of the baryon 
density of the Universe at the recombination. The Planck 
satellite (Ade et al, 2014) has produced the result 

?7cmb = (6.05 ± 0.07) x 10 -10 , (11) 


in excellent agreement with the BBN value. 

In order to compute the abundances from the BBN the¬ 
ory, the neutron lifetime r„ needs to be known with suffi¬ 
cient accuracy, because the number of neutrons available 
for fusion at Then depends on r n . According to state-of- 
the-art calculations by Coc et al. (2014), the sensitivity 
of the deuterium fraction to the values of the input pa¬ 
rameters is 


A(D/H) 

D/H 



+ 0.4 


A r n 

T n 


( 12 ) 


The parametric uncertainty due to the error on the neu¬ 
tron lifetime becomes negligible compared to the obser¬ 
vational error if AT n <C 30 s. We will explain in the next 
paragraph how the needed sub-percent accuracy has been 
achieved in the 1990’s. 


The measurement of the neutron lifetime. The 

neutron decays due to the weak interaction into a proton, 
an electron and an antineutrino with a period of about 
ten minutes. It is the simplest case of nuclear beta decay. 
Given its importance in cosmology and particle physics, 
the neutron lifetime has been measured by more than 
20 experiments. There are two distinct experimental ap¬ 
proaches to measure the neutron lifetime r n . 

1. The beam method. A detector records the decay 
products in a well-defined part of a neutron beam. 
A neutron beam is indeed radioactive due to beta 
decay and the rate of appearance in the beam of 
the decay products is 


dN _ N 
dt r n 


(13) 


This method requires (i) a determination of the 
number of neutrons N in the beam, that is, an ab¬ 
solute determination of the neutron flux and (ii) a 
detector for decay products with a well calibrated 
registration efficiency to measure dN/dt. 

2. The bottle method. A bottle is filled with UCNs. 
After a certain waiting time t , the trapped neutrons 
are emptied into a detector to count the number of 
remaining neutrons N(t). One repeats this opera¬ 
tion for various times t and the storage time r s t or 
is extracted from a fit of the storage curve 

N(t) = A(0)e" t/Tst ". (14) 


The efficiency of the neutron counter does not 
need to be known accurately, but the losses due to 








ILL, beam 1990 (corr. 1996) 



Gatchina, gravitrap 11992 

ILL, double bottle 1 1993 , u 

.Q. 


ILL. double bottle II 2000 (corr. 20121 - 

■ -O -■ 


ILL, gravitrap II 2005 - j . 



ILL. mambo II 2010 H 




865 870 875 880 885 890 895 

neutron lifetime [s] 

FIG. 1 Results of the latest neutron lifetime measurements 
with the beam method (#), material bottle method (■) and 
magnetic bottle method (<C>). Dashed lines correspond to mea¬ 
surements which have been withdrawn or corrected later on. 


absorption or inelastic scattering during collisions 
with the walls must be carefully controlled in order 
to determine r n from the storage time r stor . 

The current situation of the measurements is summa¬ 
rized in Fig. 1. References to the publications of each ex¬ 
periment can be found in the recent review by Wietfeldt 
and Greene (2011), except for the recent reevaluation of 
the latest beam result by Yue et al. (2013). 

Concerning the beam method, one can choose to mea¬ 
sure either the electron or the proton activity of the beam 
(it would be silly to try to measure the neutrino activ¬ 
ity). Although early experiments (not represented in Fig. 
1) detected the electrons coming out of the beam, mod¬ 
ern measurements count the protons. Combining the two 
most recent measurements we obtain 

r^m = (888.0 ±2.1) s. (15) 

Let us now come to the storage method where one 
measures the storage time r s t or of UCNs in a material 
bottle, which is a combination of the beta decay life¬ 
time and the lifetime due to losses at wall collisions : 
1 /7~ s tor = l/r n + 1/T wa n. Most of the time neutrons are 
specularly reflected at wall collisions, but there is a small 
probability of capture or up-scattering. The loss proba¬ 
bility can be of the order of n ~ 10~ 5 , using walls coated 
with hydrogen-free oil such as fomblin for example. In a 
typical material bottle with mean distance A « 30 cm, 
UCNs with velocity v sa 3 m/s will collide with the 
walls at a frequency / = v/X ss 10 Hz. Then we ex¬ 
pect r wa ii = 1 /ff/, ~ 10 4 s, which corresponds to a 10 % 
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correction of the beta decay lifetime. The usual strategy 
to control this effect is to measure the storage time for 
various geometries of the bottle, with different volume to 
surface ratio, and extrapolate to the ideal case of van¬ 
ishing wall collision frequency. In some experiments the 
bottle is cooled down to reduce the neutron losses due to 
up-scattering. The five UCN measurements performed 
at the ILL produced the combined value 

r^ ottle = (879.6 ± 0.8) s. (16) 

There is a 3.8er discrepancy between the beam result 
(15) and the bottle result (16) that needs to be resolved. 
To alleviate the issue of the losses at wall collisions which 
are not fully understood yet, new projects aim at confin¬ 
ing UCNs in a magnetic bottle. The neutron is a spin 
1/2 particle with a magnetic moment /j, n = —60 neV/T. 
Thus a magnetic “wall” of the order of 1 T acts as a 
repulsive potential of height 60 neV for low field seekers 
neutrons (i.e. neutrons with spin parallel to the magnetic 
field). The measurement of the neutron lifetime is still 
an active field of research. Several teams in Europe and 
in the US are currently attempting to obtain a reliable 
measurement with an accuracy of 0.1 s. 

D. The matter antimatter asymmetry and the neutron 
electric dipole moment 

Apparently our Universe is made up of matter, not 
antimatter. Not a single complex antinucleus like an¬ 
tihelium has been detected in cosmic rays. No excess 
of gamma radiation resulting from the annihilation of 
antimatter with matter is observed. Most cosmologists 
believe that the matter dominance over antimatter ex¬ 
tends to at least the whole visible Universe. In fact, the 
imbalance is tiny. When the baryons were in equilib¬ 
rium with the rest of the plasma in the early Universe 
just before the QCD phase transition, for every billion 
baryon-antibaryon pair there was one spare baryon. It 
should be noted that we quantify the matter-antimatter 
asymmetry by the baryon asymmetry. For sure there are 
also more electrons than antielectrons in the Universe to¬ 
day. However, an excess of antineutrinos over neutrinos 
in the cosmic neutrino background could perhaps com¬ 
pensate for the electron-positron asymmetry. Therefore 
we do not know for certain that an excess of matter over 
antimatter exists for leptons. As we have seen, the asym¬ 
metry parameter r) = rib /n 7 deduced from considerations 
about the primordial nucleosynthesis (10) agrees with the 
one deduced from the microwave background (11). Given 
that the two methods rely on two completely separated 
epochs in the Universe, the agreement is remarkable. 

The Sakharov conditions to generate the 
baryon asymmetry. We have solid evidence that 
the baryon asymmetry exists since before the Big Bang 


nucleosynthesis. Is it merely an initial condition of the 
Big Bang tuned to allow intelligent life to emerge? In 
the context of inflation this idea is no longer plausi¬ 
ble, because inflation would have tremendously diluted 
the initial baryon density. Well before inflation was in¬ 
vented, Sakharov (1967) proposed that the asymmetry 
could have been generated dynamically from an initially 
symmetric state. He outlined three conditions that are 
necessary for this to be possible: 

1. The baryon number should not be conserved. 

2. The Universe should at some time depart from ther¬ 
mal equilibrium. 

3. The discrete symmetries C (charge conjugation) 
and CP (charge conjugation combined with parity 
transformation) should be violated. 

See Bernreuther (2002) for a pedagogical introduction 
to baryogenesis, the hypothetical process in the early 
Universe that meets all these conditions to generate the 
baryon asymmetry. 

Electroweak baryogenesis. The second Sakharov 
condition suggests that the baryon asymmetry was gen¬ 
erated during a phase transition in the early Universe be¬ 
fore the primordial nucleosynthesis. Electroweak baryo¬ 
genesis, a mechanism that biases the baryon number 
during the electroweak phase transition, is one of the 
most attractive among the many proposed realizations of 
baryogenesis. In fact, all three Sakharov conditions are 
qualitatively fulfilled by the standard electroweak the¬ 
ory although it fails quantitatively to predict the correct 
baryon asymmetry. Nevertheless, electroweak baryogen¬ 
esis is still a viable scenario in extensions of the Standard 
Model (SM) of particle physics. New physics is required 
at or close to the electroweak scale, this makes the sce¬ 
nario testable - and falsifiable - by current or planned 
experiments. Let us take a closer look at the Sakharov 
conditions in the context of the electroweak phase tran¬ 
sition. 

Surprisingly, the SM does accommodate baryon num¬ 
ber violation, induced by non-perturbative effects asso¬ 
ciated with the nontrivial structure of the SU(2) gauge 
fields vacuum. We do not observe baryon number vi¬ 
olation in the laboratory because the process requires 
a quantum tunneling through a large energy gap with 
an extremely small probability. However, Kuzmin et al. 
(1985) discovered that baryon number violation processes 
called sphalerons were frequent in the early Universe 
when the temperature was high enough to overcome the 
energy barrier. In fact, the first Sakharov condition does 
not demand physics beyond the Standard Model. 

Next, the second Sakharov condition deserves a dis¬ 
cussion about the electroweak phase transition, when the 
Higgs field acquired a non-zero expectation value. When 




<|>/GeV 


FIG. 2 Effective Higgs potential for different temperatures 
at the electroweak phase transition. Top: calculation within 
the SM, a smooth cross-over occurs at T = 170 GeV. Bot¬ 
tom: calculation within the q !> 6 extension, a first order phase 
transition occurs at T = 120 GeV. 

coupled to the thermal bath at temperature T, the effec¬ 
tive potential of the SM Higgs field 4> can be written as 
a sum of the bare potential and a thermal correction: 

V eff (^T) = -^ 2 + ^ 4 +Vth(0,T). (17) 

At zero temperature the thermal correction to the po¬ 
tential is absent (V t h = 0) and the potential (17) has 
the famous Mexican hat shape with a minimum at </> = 
\J\x 2 /A = 246 GeV. At high temperature the thermal 
correction Vth 7^ 0 modifies the potential in such a way 
that the minimum of the potential is </> = 0. When the 
temperature decreased, the field <j> condensed from the 
symmetric state <j> = 0 to the state (f> = 246 GeV, break¬ 
ing spontaneously the electroweak symmetry. Successful 


baryogenesis requires that the phase transition be first 
order: the field must change discontinuously from <f> = 0 
to (j) > 0. Figure 2 (Top) shows the modification of the 
potential close to the phase transition at T = 170 GeV. 
The formulas for the effective potential are taken from 
the thesis of Fromme (2006), where I specified the value 
of the Higgs mass to rrih = 126 GeV conforming with the 
recent discovery at the Large Hadron Collider. We see 
that the minimum of the potential smoothly moves away 
from zero, there is no sufficient departure from thermal 
equilibrium and no baryogenesis in the Standard Model. 

Now, a non-standard Higgs potential could lead to a 
first order electroweak phase transition. As an example, 
we could add to the potential (17) a non-renormalisable 
operator of the form ^</> 6 . We plot in Fig. 2 the modifi¬ 
cation of the effective potential in the case A = 700 GeV, 
using again the formulas in Fromme (2006). In this case 
the minimum of the potential suddenly changes from 
(f> = 0 to (f> = 140 GeV at the critical temperature of 
T cr = 120 GeV, this is a first order phase transition. If 
such a transition did take place in the early Universe, 
bubbles of true vacuum started to nucleate and expand 
in a sea of false vacuum. In such a boiling environment, 
the Universe could have been driven sufficiently out of 
equilibrium for baryogenesis to be possible. 

A first order electroweak phase transition can be real¬ 
ized in more sophisticated extensions of the scalar sector 
of the SM, such as the two-Higgs-doublet model or su¬ 
persymmetric extensions. Since these models modify the 
physics at the electroweak scale, they are testable at par¬ 
ticle colliders. . 

Last, electroweak baryogenesis requires CP violation at 
the electroweak scale. It turns out that the CP violation 
contained in the SM, from phase S of the quark mix¬ 
ing matrix, is not strong enough to generate the baryon 
asymmetry. Thus CP-violating new physics is required 
to satisfy the third Sakharov condition as well. Such 
new physics is best probed by low energy precision ex¬ 
periments searching for electric dipole moments (EDM) 
of particles. 

In summary, the failure of the SM to allow for the 
electroweak baryogenesis is a hint toward the presence 
of new physics lying just above the electroweak scale 
that may be discovered both by collider experiments and 
EDM searches. 

CP violation and neutron EDM. The existence of 
a nonzero EDM for a spin 1/2 particle such as the neu¬ 
tron would imply the violation of the CP symmetry. In 
the Standard Model, the induced neutron EDM expected 
from the S phase is tiny, d n « 1CU 32 e cm. This value is 
to be compared to the current best limit obtained at the 
ILL by Baker et al. (2006): 

\d n \ < 3 x 10 -26 ecm (90% C.L.). (18) 

Therefore, improvements of the neutron EDM measure- 
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FIG. 3 One loop diagram contributing to the quark EDM. 
The quark EDMs then transfer almost directly to the neutron 
EDM. 


ment are motivated by the potential discovery of a new 
source of CP violation beyond the SM. As for any low en¬ 
ergy observable, the new physics at a high energy scale 
manifests itself via quantum loops involving virtual par¬ 
ticles. Figure 3 shows a Feynman diagram contributing 
to a quark EDM via the CP-violating vertices of heavy 
scalars and fermions with masses M. 

Generically the neutron EDM induced by such a loop 
amounts to (Pospelov and Ritz, 2005) 

d n ss ( ^ J x sin (j) x 10~ 25 ecm (19) 

where M is the mass of the particles running in the loop. 
In this case the heavy particles couple strongly with the 
quark (such as the SUSY coupling between quark, squark 
and gluinos) with a CP-odd vertex multiplied by sin </>. 
The CP-odd part usually originates from the imaginary 
part of some parameter in the Lagrangian and (j) would 
then correspond to the CP-violating phase of that spe¬ 
cific parameter. Thus, considering natural CP violation 
(sin ip « 1) in the new heavy sector, the neutron EDM is 
sensitive to new physics at the multi-TeV scale. 

Searches for permanent electric dipole moments of 
other particles (protons, deuterons, muons, atoms, 
molecules) are complementary probes of CP violation 
above the electroweak scale. Many experimental efforts 
are under way, they have been compiled recently by Kirch 
(2013). Below we concentrate on the search for the neu¬ 
tron EDM. 

Measuring the neutron EDM. The neutron EDM 
measurement is based on the analysis of the Larmor pre¬ 
cession frequency of neutrons, stored in a volume perme¬ 
ated with electric and magnetic static fields either paral¬ 
lel or antiparallel. For such configurations, the precession 
frequency f n reads 


hf n = n n B ± d n E (20) 

where and d n are the magnetic and electric dipole 
moments of the neutron and h is Planck’s constant. The 



FIG. 4 The nEDM apparatus installed at the PSI. The pre¬ 
cession chamber sits in a large cylindrical vacuum chamber 
(with front end open on the picture), surrounded by a four- 
layer mumetal magnetic shield. The UCN guide injects the 
neutrons from the bottom. High voltage is applied using the 
HV feed-trough connected to the top electrode. 


frequency difference of these two configurations gives 
directly access to the neutron EDM: d n = hAf n /AE. 
To measure the precession frequency, we use Ramsey’s 
method of separated oscillatory fields which provides a 
statistical precision per cycle of 


cr(d n ) 


h 

2 aETVN 


( 21 ) 


where T is the precession time, a is the visibility (re¬ 
lated to the polarization of the neutrons) and N the to¬ 
tal number of detected neutrons. Using stored ultracold 
neutrons, the precession time can be set to T = 200 s, 
five orders of magnitude longer as compared to the first 
experiment performed by Smith et al. (1957) using a neu¬ 
tron beam! 

The main experimental challenge in current experi¬ 
ments consists in achieving a magnetic field homogeneity 
at the level of 10 -4 over a volume of typically 20 i while 
maintaining a temporal stability of about 10~' over 100 
s. These requirements are necessary to control the sub¬ 
tle systematic effects (see for example Pignol and Roccia 
(2012)). Atomic magnetometry and magnetic shielding 
techniques are therefore at the core of such a measure¬ 
ment. 

As of 2015 there are six projects worldwide competing 
for improving the measurement of the neutron EDM, all 
at different stages. 


1. A PNPI group is currently operating a room- 
temperature double-chamber spectrometer at the 
ILL (Serebrov et al, 2014). Later the apparatus 
will be moved back to PNPI where a new UCN 
source will be built. 
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2. Another PNPI group pursues a completely differ¬ 
ent experimental method. It exploits inter-atomic 
electric fields in crystals that are 10 4 times stronger 
than electric fields produced by charging electrodes 
as in all other experiments. However the interac¬ 
tion time is much shorter because it is a neutron 
scattering technique. The result obtained recently 
at the ILL (Fedorov et ai, 2010) is two orders of 
magnitude less precise than the conventional UCN 
method. 

3. A German project at Munich (Altarev et ai, 2012) 
is waiting for the construction of a new UCN source 
at the reactor FRM-2. The newly built magnet¬ 
ically shielded room performs very well (Altarev 
et al . , 2014). This large room is designed to host a 
double-chamber room-temperature experiment. 

4. A Canada-Japan collaboration plans to build a su- 
perthermal UCN source and a room-temperature 
EDM apparatus at TRIUMF Vancouver, based on 
the prototype operated at RNCP Osaka (Masuda 
et al., 2012, 2014; Matsuta et ai, 2013). 

5. An ambitious project is pursued by the US commu¬ 
nity based on the original proposal of Golub and 
Lamoreaux (1994). The idea is to run the exper¬ 
iment in superfluid 4 He, that serves both as a su- 
perthermal UCN source and a precession chamber. 
Traces of spin-polarized 3 He will be injected in the 
volume to both detect the neutrons and act as a co¬ 
magnetometer (Ye et al, 2009). The project is still 
in the R&D phase and will start real data taking 
at SNS Oak Ridge in 2020 at best. 

6. A European collaboration (Belgium, France, Ger¬ 
many, Poland, Switzerland and UK) is running an 
experiment at the PSI UCN source (Baker et al . , 
2011). We are currently operating an upgraded ver¬ 
sion of the RAL-Sussex spectrometer (see picture 
Fig. 4) which has produced the lowest experimen¬ 
tal limit for the neutron EDM and that we moved 
from the ILL to the PSI in 2009. In the same 
time we have started the design of n2EDM, a next 
generation room-temperature double-chamber ap¬ 
paratus in view of its delivery around 2018. With 
the current apparatus we aim at slightly improv¬ 
ing the present limit on the neutron EDM whereas 
with n2EDM we aim at a precision of better than 
1(T 27 e cm. In addition, we have produced sev¬ 
eral spin-off scientific results with the apparatus 
(i) a search for neutron to mirror-neutron oscilla¬ 
tions (Altarev et al, 2009b; Ban et al, 2007) (ii) a 
sensitive test of Lorentz invariance (Altarev et al., 
2009a, 2010, 2011) (iii) a measurement of the neu¬ 
tron magnetic moment with an uncertainty of 0.8 
ppm (Afach et al . , 2014) (iv) a search for Axionlike 
particles (Afach et ai, 2015). 


Although the experiments are getting more and more 
difficult (in the 1950s and 1960s two or three persons 
could run an experiment, nowadays author lists rarely 
count less than 20 people) prospects are good to improve 
the accuracy on the neutron EDM by more than an order 
of magnitude in the next decade. 


In conclusion, experiments with neutrons address the 
cosmological question of the origin of the baryonic mat¬ 
ter. Today we understand how protons and neutrons 
combined in the early Universe to form nuclei. Experi¬ 
ments measuring the neutron lifetime contributed a fare 
share of the successful prediction of the Big Bang nucle¬ 
osynthesis. However we still do not know how the neu¬ 
trons and protons were created in the first place, because 
we do not understand the baryon asymmetry of the Uni¬ 
verse. Planned measurements of the neutron EDM will 
contribute to solve this question. 
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II. CHAMELEON DARK ENERGY 

The 2011 Nobel Prize in physics was awarded to S. 
Pelmutter, B. Schmidt and A. Riess for the discovery of 
the late acceleration of the expansion of the Universe. 
Although the acceleration is now established as a fact, it 
remains a great puzzle. We do not know the nature of 
the Dark Energy responsible for the acceleration. 

The two other cosmological puzzles, the origin of the 
baryon asymmetry and the nature of the Dark Matter, 
are thought to be related to physics at high energy, be¬ 
yond the electroweak scale. On the contrary, the en¬ 
ergy scale associated to the Dark Energy has the peculiar 
value of 2 meV. It suggests that there is still new physics 
in the infrared to be discovered, perhaps with low energy 
precision experiments using neutrons? 


A. The accelerated expansion of the Universe 

A linear relation between the distance and radial ve¬ 
locity among galaxies, v = H 0 d , was first established 
by Hubble (1929). The first Hubble diagram is shown 
in Fig. 5. Incidentally, due to an error in the dis¬ 
tance evaluations, Hubble derived an expansion rate of 
500 km/s/Mpc which is much larger than the most recent 
determination by Ade et al. (2014): 

H 0 = 67.3 ±1.2 km/s/Mpc. (22) 

Hubble’s linear law was discovered using observations 
of nearby galaxies. At much larger distances, the rela¬ 
tion is not expected to be linear anymore. The expan¬ 
sion of the Universe was thought to slow down because of 
the attractive power of gravity. The deceleration of the 
Universe was then searched for by drawing a Hubble di¬ 
agram with more distant objects. Explosions of type IA 
supernovae are (i) very bright, they can be seen at cos¬ 
mological distances (ii) very reproducible, they all have 
almost the same intrinsic brightness. A sufficient number 
of supernovae observations using the Hubble space tele¬ 
scope and ground based telescopes were reported by two 
teams (Perlmutter et al, 1999; Riess et al., 1998). The 
supernovae Hubble diagram of one team is shown in Fig. 
5. 

Surprisingly, it was found that the expansion of the 
Universe accelerates. Let us analyse this incredible find¬ 
ing in the framework of a flat Universe whose expansion 
is driven by an homogeneous substance with an energy 
density p and a pressure p (i.e. a perfect fluid). We as¬ 
sume an equation of state of the substance of the form 
p = tup. From the Friedmann equations (4), (5) one de¬ 
rives the acceleration: 

d 47 tG 

- =-—p(l + 3w). (23) 

fl O 

A substance made up of non-relativistic matter (w = 0) 
or radiation (w = 1/3) would decelerate the expansion. 



Velocity-Distance Relation among Extra-Galactic Nebulae. 





FIG. 5 Top: diagram taken from Hubble (1929) representing 
the velocity (deduced from the redshift) of nearby galaxies as 
a function of distance. Bottom: supernovae Hubble diagram 
taken from Perlmutter et al. (1999) representing the magni¬ 
tude (indicating the distance) as a function of the redshift. 
Notice that the highest velocity of 1000 km/s in the 1929 di¬ 
agram corresponds to the redshift z = 0.003 which is smaller 
than the smallest redshift in the supernovae diagram. 


Instead, a positive acceleration requires w < —1/3. The 
mysterious substance driving the expansion must have a 
negative pressure. This substance is generically called 
Dark Energy. 


B. The ACDM concordance model 

In the standard model of cosmology, the acceleration 
of the expansion is attributed to a perfect fluid with w = 
— 1. According to the second Friedmann equation (5), 
the energy density of this fluid is constant. It does not 
dilute with the expansion, it is a cosmological constant. 

In the late Universe, the expansion is driven by two dis¬ 
tinct perfect fluids: the cosmological constant and some 
non-relativistic matter (radiation became unimportant a 
few million years after the Big Bang). Hence the name 
“ACDM” for the standard model, which stands for a cos- 
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FIG. 6 The age of the Universe calculated in the flat ACDM 
model as a function of Da- The vertical orange band repre¬ 
sents the Dark Energy fraction reported by Ade et al. (2014). 


mology with a cosmological constant and cold dark mat¬ 
ter (baryonic matter accounts for only a fraction of the 
non-relativistic matter). 

We note pm,o and Pa,o the present value of the energy 
density associated with the matter and the cosmological 
constant, and p Ci o = the critical density. We use the 
standard notation !1 m = Pm,o/Pc,o and Da = Pa,o/Pc,o> 
which represent the fraction of the total energy density in 
the form of matter and cosmological constant. For a flat 
Universe, Eq. (4) considered at the present epoch implies 
Hm + Da = 1. From the second Friedmann equation (5) 
the past values of the densities are 

PAf(t) = , PA{t) = pa,o- (24) 

a(t) 3 

Then, Eq. (4) can be brought to the form 

(i) 2 = "»(^ + nA )- (25) 

This is the equation describing the expansion rate in the 
flat ACDM model. Let us calculate the age of the Uni¬ 
verse, by transforming the previous equation into: 


The age of the Universe is plotted in Fig. 6 as a func¬ 
tion of Da- We can compare it to the age of the oldest 
stars, estimated to at least 11 Gyr (Bartelmann, 2010). 
The Einstein-de Sitter model, i.e. the flat matter dom¬ 
inated Universe with Dm = 1 and Da = 0, is excluded 
because the Universe cannot be younger than the stars. 
The long-standing disagreement between the Hubble ex¬ 
pansion rate and the age of the oldest stars has finally 
been resolved in the ACDM model with the introduction 
of the cosmological constant. 

In addition to the supernovae and the age of oldest 
stars, the acceleration of the Universe is now supported 
by several other independent cosmological observations 
including baryon acoustic oscillations, weak gravitational 
lensing and counts of galaxy clusters. These observa¬ 
tional probes of cosmic acceleration were recently re¬ 
viewed in Weinberg et al. (2013) and future prospects 
are described. 

Since the discovery of the accelerated expansion, the 
precision of the cosmological data increased spectacu¬ 
larly, in particular the measurements of the CMB. The 
ACDM model, which emerged at the end of the last cen¬ 
tury as a concordance cosmology, is still able to describe 
all observations. ACDM fits indicate that the energy 
budget of the Universe today is: 

Dft = 5 % D dm = 27 % D A = 68%, (28) 

where the non-relativistic matter has a baryonic and non- 
baryonic component (Dm = Ddm + Db). These three 
components driving the expansion of the Universe today 
are associated with the three major puzzles in cosmology. 
It is often said that the baryonic content of the Universe 
is the only part which is well understood. In fact, the 
very presence of baryons results from the initial matter- 
antimatter asymmetry of the Universe, it cannot be ex¬ 
plained by the standard theory. Next, the microscopic 
description of the non-baryonic Dark Matter is lacking. 
It could be made up of massive weakly interacting parti¬ 
cles a.k.a. WIMPS. It could also well be something com¬ 
pletely different like a Bose-Einstein condensate of very 
light scalar particles like Axions. Last, the microscopic 
nature of the Dark Energy is completely unknown. Con¬ 
trary to Dark Matter, no well motivated extensions of 
the Standard Model of particle physics provide a natural 
candidate for the Dark Energy. 


a y/a 

+ « 3 Da 


(26) C. The cosmological constant problem 


Next, we integrate the equation from t = 0 to t = to (to is 
the age of the Universe). Recall that in our conventions 
a(t 0 ) = 1. We get 

1 /' 1 ^fada 2 asinh^/DA/DM 

° H 0 J 0 v^m + o 3 Da 3if 0 VQa 

(27) 


In natural units (h = c = 1), the energy density asso¬ 
ciated with the Dark Energy is 

PA = Da/ 3 c ,o = (2.2 meV) 4 . (29) 

The fact that the energy scale associated with the Dark 
Energy is so curious constitutes the cosmological con¬ 
stant problem. In the ACDM model, the Dark Energy is 
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simply the energy density of the vacuum, or equivalently 
a cosmological constant. In principle the vacuum energy 
is expected to receive contributions from several sources. 

At the classical level, any scalar field </> contributes to 
the vacuum energy by the amount V(4> o) where V(</>) is 
the energy density function of the field and </>o is the vac¬ 
uum value of the field. One could argue that the zero en¬ 
ergy is arbitrary, we can choose to define it as the state of 
lowest energy for all fields in the Universe. All local phe¬ 
nomena do not depend on this choice of the zero energy, 
there is no reason that it should have such a dramatic 
gravitational effect on cosmological scales. However, in 
the case of phase transitions, the effect of zero energy 
cannot be ignored by simply declaring that the zero en¬ 
ergy cannot gravitate. During the electroweak symmetry 
breaking mentioned in section I.D, the energy density 
associated with the Higgs field changes, by an amount 
as large as (100 GeV) 4 . The vacuum energy was much 
larger before the electroweak phase transition. This dif¬ 
ference in vacuum energy is definitely a physical mean¬ 
ingful quantity, according to straightforward application 
of the general relativity theory in a purely classical con¬ 
text. If differences in vacuum energy certainly gravitate, 
what determines the actual value of the vacuum energy? 
It could perfectly correspond to an intrinsic property of 
the Universe that should be treated as another free pa¬ 
rameter. Why then is it adjusted so that the vacuum 
energy after the electroweak phase transition is almost 
zero, but not quite? It seems that Nature decided to 
“lift” the Mexican hat potential of the Higgs field in a 
very precise, fine-tuned amount. 

In addition, the vacuum energy receives contributions 
from quantum fluctuations induced by virtual particles. 
Quantum field theories predict that every bosonic de¬ 
gree of freedom of frequency w has a zero-point energy 
of Hu>/2. Since there are in principle an infinite number 
of degrees of freedom in the quantum fields, the total 
vacuum energy is technically infinite, unless the contri¬ 
bution from the fermionic and bosonic degrees of freedom 
compensate each other. If physics can be described by 
an effective local quantum field theory up to the Planck 
scale, 


m P i = ^— = 2.4 x 10 18 GeV, (30) 

dimensional analysis indicate that the vacuum energy in¬ 
duced by quantum fluctuations should be of the order of 
TOpj ~ 10 12 °pa- This is definitely the worse theoretical 
prediction ever. Supersymmetric compensation between 
the fermionic and bosonic quantum fluctuations above 
the SUSY breaking scale could reduce the disagreement 
between the expected and observed vacuum energy from 
10 120 down to 10 58 , still enormous. 

In summary, the vacuum energy can be separated as 
a sum of two terms: pa = Pa, class + PA,quant- The term 


induced by quantum fluctuations pa, quant cannot be cal¬ 
culated in a self consistent way in quantum field theory 
but it is expected to be huge from dimensional analysis. 
The classical term pa, class is the energy density of the 
scalar fields in the ground state, it is a free parameter of 
the theory. At the end, the vacuum energy can perfectly 
be considered as a free parameter in the present incom¬ 
plete theories of particle physics and gravity. However, 
there is a naturalness problem very similar to the hierar¬ 
chy problem of particle physics in the sense that “bare” 
value has to be extraordinarily fine-tuned to almost com¬ 
pensate for the huge quantum contribution. 

Well before the discovery of the accelerating Universe, 
Weinberg (1987) discussed an anthropic explanation to 
the cosmological constant problem. He gave an upper 
bound on the vacuum energy density: 

PA < (10meV) 4 . (31) 

The argument goes as follows. If the vacuum energy is 
too large, the accelerated phase of the expansion starts 
too early. Then the inhomogeneities of the mass den¬ 
sity do not have time to grow sufficiently before getting 
stretched away by the accelerated expansion. Only if the 
Weinberg bound (31) is satisfied can the perturbations 
grow to form stars and galaxies. Otherwise the Universe 
would consist in an homogeneous fluid in accelerated ex¬ 
pansion, quite unfriendly for life. Imagining that many 
universes exist, with a different value for the cosmologi¬ 
cal constant in each, it is not surprising that we happen 
to live in a seemingly very special universe where life is 
possible at all. Indeed the observed vacuum energy den¬ 
sity (29) is not particularly fine-tuned with respect to the 
range of possible values in the anthropic sense given by 

(31) U 

This anthropic explanation is further supported by the 
landscape of string theory: there is a huge number of pos¬ 
sible vacua of the string theory associated with the many 
many different ways to conrpactify the extra dimensions. 
One might believe that all these versions of string the¬ 
ory “exist” somehow, and intelligent life can emerge only 
in a tiny fraction of a vast number of possible versions. 
The fine-tuning of the parameters could then be just an 
illusion. 

In addition, the eternal inflation scenario supports the 
plausibility that many worlds can “exist” in a multiverse. 
In these theories, the Universe is inflating due to the en¬ 
ergy density of a scalar field - the inflaton - which is in 
a false vacuum. The inflaton decays in the true vacuum, 
stopping the inflationary phase. But this process hap¬ 
pens only locally, forming a bubble of true vacuum sur¬ 
rounded by the rest of the Universe which is still inflating. 
In some versions of the theory, the part of the Universe 
which is still inflating grows fast enough to accommo¬ 
date for the continuous formation of bubbles (the speed 
of the bubble walls could be slower than the exponential 
expansion of the false vacuum space). In this case the 
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inflationary phase of the Universe lasts forever in some 
regions of the Universe and causally disconnected bub¬ 
bles are produced on and on. Different bubbles having 
different cosmological constants is not strictly speaking 
a prediction of eternal inflation, but if inflation is due to 
Planck-scale physics or stringy effects, who knows? We 
might be living in one out of a vast number of bubbles, 
where the vacuum energy density is such that it does not 
forbid the appearance of intelligent life. 


D. The quintessence as Dark Energy 

The anthropic explanation of the cosmological con¬ 
stant problem is both elegant and fascinating. However 
it rests on speculative ideas, to say the least. This ex¬ 
planation is a source of great inspiration for imagination 
and philosophy, unfortunately it gives poor insight for 
further observations and experiments. It is thus neces¬ 
sary to explore other possibilities. We can take the view 
that Dark Energy is a dynamical process revealing new 
physics at the energy scale p « 1 meV. A review of the 
possible dynamics of Dark Energy and its observational 
consequences can be found in Copeland et al. (2006) . Al¬ 
though none of these models really solve the cosmologi¬ 
cal constant problem in a top-down way, they are useful 
bottom-up approaches to guide future observations and 
experiments to refine our knowledge of Dark Energy. 

The most popular routes to a dynamical extension 
of the simple ACDM parametrization are quintessence 
models. Akin to inflationary models, the accelerated ex¬ 
pansion is due to the energy density of a scalar field. 
The scalar field responsible for the acceleration of the 
Universe during the inflation era is called the infla- 
ton, whereas the hypothetical scalar field making up the 
Dark Energy driving the late acceleration is called the 
quintessence. There are even models explaining both the 
early inflation and the late acceleration of the Universe 
using a single scalar field (Peebles and Vilenkin, 1999). 

Let us review the theory of the gravitational dynam¬ 
ics of a scalar field. The action for gravity and the 
quintessence field ip with the so-called minimal coupling 
is 

S = J d^xV^g 7 l + (32) 

where g = det(g^) is the determinant of the metric ten¬ 
sor, 1Z = g^ v Rfj, v is the trace of the Ricci curvature tensor 
Rfiv and 

£<p = - V{y) (33) 

is the Lagrangian of the scalar field with a generic poten¬ 
tial V{ip). 

Einstein’s equations governing the response of the met¬ 
ric to the scalar field are obtained from the variational 


principle SS/6g flLI = 0: 


Gpv = 8 ttGT^ v , (34) 


where = R^ 


rptp _ O 

Sg^ 


l/2g liu n is Einstein’s tensor and 
- g^vC-v (35) 


= d^ipdvip 


1 

2 


g^g a0 d a (pdg(p + g^V(ip) 


(36) 


is the stress-energy tensor of the scalar field. 

In the homogeneous Universe, the scalar field is uni¬ 
form in a Robertson-Walker metric (3). The problem 
reduces to only two degrees of freedom a(t) and <p(t) and 
tensors take a diagonal form: 


9nu = diag (l,-a 2 ) (37) 

= diag (1,-1/a 2 ) (38) 

= diag (sH 2 ,a 2 [-3H 2 - 2 H]\ (39) 

T £v = diag Qy > 2 + V{tp),a 2 \^g> 2 - V(y>)]^ (40) 


The 00 component of Einstein’s equations reduces then 
to the Friedmann equation 


h2 = 8 ?(^ 2+vw )' (41) 


In addition, the evolution of the scalar field is given by 
the variational principle 5S/5tp = 0, which in the homo¬ 
geneous Universe reduces to 


<p + 3H<p + V(<p) = 0. (42) 


The coupled dynamics of the expansion (i.e. the time 
evolution of H = a/a and tp) is completely specified by 
the two equations (41) and (42). It is useful to notice that 
these two equations are equivalent to the two Friedman 
equations (4), (5) with p = ip 2 / 2+V(v?) the energy density 
of the quintessence andp = <p 2 /2—V(<p) its pressure. The 
equation of state parameter 


P = y> 2 /2 - V{ip) 
p <p 2 /2 + V(ip) 


(43) 


is a dynamical quantity that can evolve with time. 

Different cosmologies can be built depending on the 
potential V (</?). If <p corresponds to a standard mas¬ 
sive scalar field, such as the Axion, the potential V(<p) = 
\p 2 p 2 is parabolic. For the harmonic potential there is 
equipartition between the kinetic energy p 2 /2 and the 
potential energy Vipp) and therefore w = 0. In this case 
ip is a candidate for pressurless cold Dark Matter. It is 
also possible to build models which mimic a cosmological 
constant with w « — 1, where the field does not develop 
significant kinetic energy. Such a model has been first 
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proposed by Ratra and Peebles (1988), with the poten¬ 
tial given by 


V( l p) = Mi(^y (44) 

where M\ is a new mass scale and n is a positive exponent 
called the Ratra-Peebles index. 

As a generic feature of quintessence models, a de¬ 
parture from w = —1 is expected. Therefore there is 
a chance that future supernovae surveys could detect 
a deviation from the ACDM model where w = — 1. 
Now, is there another way to reveal the presence of the 
quintessence field <p? If the field were coupled to ordinary 
matter, then it would mediate a fifth force. It is then very 
appealing to design laboratory experiments searching for 
this new force. 


E. The chameleon 

A quintessence field coupled with matter mediates a 
new force with a range extending up to cosmological 
scales. At first it was generally assumed that the coupling 
must be very small, much smaller than the gravitational 
strength, in order to satisfy the stringent tests of the 
equivalence principle, as well as the tests of general rel¬ 
ativity in the Solar system. However, Khoury and Welt- 
man (2004a,b) discovered that the coupled quintessence 
scalar theory automatically features a very efficient mech¬ 
anism to suppress the force called the chameleon mecha¬ 
nism. Then, Brax et al. (2004) explored the cosmological 
consequences of the chameleon theory to find that it is a 
nice viable Dark Energy candidate. 

Let us explain the chameleon screening mechanism. 
We start by considering the Lagrange density of the 
scalar field ip coupled with a fermion ijj (with a mass 
to) in a Minkowski spacetime: 3 

1 777 _ 

£ = -df.ip d^ip - V(ip) - /?—- ip ipi’. (45) 

2 rapi 

The dimensionless constant (3 corresponds to the strength 
of the coupling relative to gravity. From the Lagrange 
density we can deduce in principle how p affects a parti¬ 
cle, and conversely, how the presence of particles gener¬ 
ates a field. 


3 The term f}m/mp\(p originates from a conformal coupling 
of the scalar field to the fermion. The conformal coupling con¬ 
sists in replacing the metric g^ u in the fermion Lagrange density 
by a function of the scalar field A 2 (ip)g^ u : £ m {A 2 (<p)g {AW 5 VO- 
A popular choice of the conformal coupling function is A(<p) = 
e /3v?/m P1 Since the excursions of the field always satisfies 
/3ip/mp\ < 1 in concrete cosmological or laboratory situations, 
we can approximate A(ip) = 1 + (3ip/mp\. It finally produces the 
Yukawa coupling term (3m/mp\(p pif. 



FIG. 7 The black dashed line represents the runaway po¬ 
tential V (</?). The solid lines represent the effective potential 
Veff(^) where the environment density p is low (blue) and 
where the density is high (red). 


1. A fermion evolving in a field p(r) is affected by the 
potential 


m 


V(r) = /3 - ip( r), 

TOpi 


(46) 


corresponding to a force: 


TO 


F = — (3 - V<p(r). 

TOpi 


(47) 


2. The field <p(r) generated by a distribution of 
fermions is given by the Euler-Lagrange equation 


( dc \ dc = 

AdidMJ dip 


(48) 


More specifically, the Euler-Lagrange equation deduced 
from the Lagrange density (45) reads 

0 M 0'V + V , ( V >)+0— =0, (49) 

TOpi 

where p is the mass density of the fermions (we have 
assumed the non-relativistic limit mipil) = p). 

It is then apparent that, within an environment of den¬ 
sity p, the dynamics of the chameleon is not governed by 
V(<p), but rather by the effective potential 

V off =V(<p) + /3-^- ip, (50) 

TOpi 

in the sense that Eq. (49) can be written in the form 

d^ip + V' ff (<p) = 0. 

Let us recall that the equation governing a basic mas¬ 
sive scalar field (with a potential \p?ip 2 ) is the Klein- 
Gordon equation d^d^ip+p^ip = 0, where p is the mass of 







16 


the elementary excitations of the field. In this well-known 
case, the field generated by a static point source (if the 
field is coupled to fermions, this corresponds to a source 
term in the right hand side of the Klein-Gordon equation) 
takes the form of a Yukawa potential ip{r) ~ e~^ r /r. 
It means physically that the field mediates an interac¬ 
tion between fermions with a range l//i (in S.I. units, 
A = ) • A very massive field will mediate a short-range 
force, whereas a light field can travel long distances. 

We now come back to the more complicated situation 
of the coupled quintessence field. We will discuss the 
consequences of the specific features of the effective po¬ 
tential depicted in Fig. 7. Assuming that the density 
p is uniform everywhere, then there is a static and uni¬ 
form solution to Eq. (49), given by the minimum of the 
effective potential V' ff (</J m i n ) = 0. In the case of the 
Ratra-Peebles potential (44) the minimum of the effec¬ 
tive potential reads 

= (51) 

For small perturbations of the field around this uniform 
value (induced by an extra mass on top of the environ¬ 
ment density p for example), we can approximate the 
effective potential around the minimum by 

Peff(^) ~ V e ff(^min) 4“ T^min) E e ff (^min) ■ (52) 

We can attribute an effective mass p for the perturba¬ 
tions of the field, with 

p 2 = V l ' ff (<p niin ) = P——. (53) 

W p l tPmin 

The effective mass is associated with the curvature of 
the effective potential at the minimum. Clearly, what is 
happening here is that the mass of the field is density- 
dependent, in such a way that the range of the force me¬ 
diated by the field shrinks in a high density environment. 
For this reason, the scalar field is called “chameleon”: it 
adapts its properties to the environment to evade obser¬ 
vation. 

Consider for example the Sun as a source of the 
chameleon field. Since the density in the interior of the 
Sun is quite high, the chameleon is effectively massive 
in the Sun. Therefore, a test particle situated outside 
the Sun is not affected by the inner part of the Sun be¬ 
cause the field created by the mass inside the Sun cannot 
propagate to the outside. In fact, it is shown in Khoury 
and Weltman (2004a,b) that when the coupling /3 is large 
enough, then only a thin shell at the surface of a large 
body such as the Sun or the Earth contributes to the 
force on a test particle outside the body. Then, in a more 
elaborate analysis, Mota and Shaw (2007) concluded that 
strongly coupled chameleons, i.e. /3 1, are not ruled 

out by terrestrial and solar system tests of gravity. 



FIG. 8 Sketch of the one-dimensional problem. 


F. ID solutions of the chameleon equation 


In order to gain more insight about the properties 
of the chameleon field, we shall now solve explicitly 
the chameleon equation in the case of the simple one- 
dimensional problem illustrated on Fig. 8. We consider 
a plate of material (to be specific, a plate of aluminium 
with thickness 2d = 1 mm) surrounded by perfect vac¬ 
uum. 

We consider a static solution in one dimension, there¬ 
fore d^d^ip = — A <p = — d ■ In this case the 

chameleon equation (49) becomes 


d 2 tp 

dx 2 




(54) 


In addition, the symmetry of the problem implies 
p(—x) = <p{x). We will therefore solve for x > 0 with 
the boundary condition i p' 0 = 0. For simplicity we treat 
the case with n = 2 but other Ratra-Peebles indices can 
be treated similarly. 


Solution in vacuum. In the vacuum (X > 0), the 
chameleon equation reads 


fv = _ 2 < 

dX 2 ip 3 


(55) 


we can find a family of solutions with different values ips 
at the boundary: 

<p(X) = ysyJl + X/Xs (56) 


with Xs satisfying 


V^2 ip 2 s 

4 M 3 K 


(57) 


Notice the remarkable expression of the derivative at the 
boundary 


ip’s — y/2 


Ml 

Vs 


(58) 
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Next, we deal with the inside of the plate (0 < x < 
d). No exact solution exists for the problem, we need to 
consider two asymptotic regimes. 

Solution inside the plate: linear regime. We 

start by considering what we call the linear regime, where 
po ^ yViin ■ We will derive a posteriori the corresponding 
condition for the parameter /3. We recall that p m ; n is the 
value of the field minimizing the effective potential inside 
the plate, it is given by Eq. (51). In this regime, the 
effective potential is dominated by the coupling term and 
we can neglect the self-interaction term V(p). Therefore, 
the chameleon equation reduces to a Poisson equation : 

*V"f,JL = ™L. ( 59 ) 

dx TOpi p d min 

with the solution 

A/f6 

ip{x) = p 0 + - 3 ^- x 2 (60) 

r min 

We now need to connect this solution to the vacuum so¬ 
lution at the boundary by equalizing the value of the field 
and the derivative. We get the solution 

^- = —1 - (‘^P) 3 -W) 2 ( — V. (61) 

Ma V2{M A d)\M A J 

Now we can work out a posteriori the domain of validity 
of the linear approximation. We see from the solution 
that the condition po ^$> p m i n is equivalent to the condi¬ 
tion M A d < {<Pmin/M A ) 2 , or, expressed in terms of the 
basic parameters: 


the field inside the plate is very close to p m \n- We re¬ 
call that (p m i n would be the uniform value of the field 
in the situation where the plate with density p fills the 
entire space uniformly. In the saturated regime we can 
write p(x ) = ^ m in + <t>, with </> <C p m \n- In this case, the 
chameleon equation reduces to the Klein-Gordon equa¬ 
tion 


d 2 4> 
dx 2 


= P<t> 


(65) 


where the effective mass p satisfies Eq. (53), that is: 


..2 _ Q P 3 6M a 

r — P 4 

m pl ^min V^min 


( 66 ) 


At this point it is useful to notice that the opposite of 
the linear regime condition Eq. (62) can be written in 
the form: 


saturated regime : pd 1 . 


(67) 


Physically, it means that in the saturated regime the 
range of the force mediated by the chameleon is much 
shorter than the thickness of the plate. 

The solution of the Klein-Gordon equation is: 


<f>{x) = cosh(px). ( 68 ) 

The continuity of p and p' at the boundary provides a 
quadratic equation for fo: 

fop sinh(pd) ( p min + </> 0 cosh (pd)) = V2 M A , (69) 


linear regime : 



1 

(. M A d ) 3 / 2 ' 


(62) 


or equivalently, we get a quadratic equation for ps = 
Thnin + <t> 0 cosh pd: 


Let us evaluate numerically this condition in the case 
of our aluminium plate. Equations are expressed in nat¬ 
ural units, so we will express all masses in eV and all 
distances in eV -1 . We will consider the energy scale of 
the chameleon potential to correspond to the Dark En¬ 
ergy scale: 

M a = Pa /4 = 2.4 x 1CT 3 eV. (63) 

The Planck mass amounts to m v \ = 2.4 x 10 27 eV. The 
half-thickness of the plate is d = 0.5 mm = 2.5 x 10 3 eV ^ 1 
and the density of aluminium is p = 2.7 g/cm 3 = 1.2 x 
10 19 eV 4 . We find that the solution is described by the 
linear regime when 


P < 0.2. (64) 

The linear regime does not describe strongly coupled 
chameleons (i.e. /3 1 ) in the presence of ordinary ob¬ 

jects of normal size! 

Solution inside the plate: saturated regime. In 

the opposite regime, which we call the saturated regime, 


(PS - ‘Pmin)PS = V2 


cosh^id 
p sinh pd 


(70) 


Interpretation of the solution. The solution of 
the chameleon equation on the right side of the plate is 
plotted on Fig. 9, for three different values of the cou¬ 
pling constant j3. The field inside the plate is “attracted” 
to small values. Outside the plate, the field wants to grow 
to minimize the potential V(p). 

For strongly coupled chameleons, the field inside the 
plate takes a nearly uniform value p = p m i n . When 
the field is saturated inside the bulk, the plate acts as a 
screen: if a source mass moves on the left side of the plate, 
the field will be modified on the left side but it will remain 
saturated to p mi n inside the plate. Therefore the solution 
on the right side of the plate will remain unaffected. The 
plate shields the right side from what happens on the left 
side. This screening mechanism makes the chameleon 
field difficult to detect by laboratory experiments. 

Another important feature of the chameleon is the sat¬ 
uration of the field. We see in Fig. 9 that the solution 
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FIG. 9 Solution of the ID chameleon equation with n = 2. 
The half-thickness of the aluminium plate is d = 0.5 mm. 
The blue curve (0 = 0.05) is calculated using the linear ap¬ 
proximation. The black (0 = 1) and red (0 = 20) curves are 
calculated using the saturated approximation. 



FIG. 10 Gradient of the chameleon field at X = 10 pm from 
the surface of the aluminium plate, as a function of 0. The 
blue curve is calculated using the linear approximation, the 
red curve is calculated using the saturated approximation. 


outside the plate becomes independent of 0 for large val¬ 
ues of 0. To better understand this, let us discuss the 
gradient of the field at the immediate vicinity of the 
plate, say X — 10 pm. The gradient of the field is anal¬ 
ogous to an electric field (which is the gradient of the 
electric potential). The force acting on a test mass is 
proportional to in the same way as the force acting 
on a charged particle is proportional to the electric field. 
There is no saturation effect for the electric field, in the 
sense that an electrically charged plate produces at the 
vicinity of the surface an electric field proportional to the 
charge density of the plate. Fig. 10 shows the gradient 
of the chameleon field as a function of 0. In the linear 
regime, the gradient grows linearly with 0 (or p) as for 
the electric field (hence the name of the linear regime). 
On the contrary, in the saturated regime, the gradient 
saturates to a finite asymptotic value. 

As a matter of fact, when searching for strongly cou¬ 
pled chameleons, one should think of solid matter as a 
zero boundary condition for the field ip. In this case the 
field in the vicinity of the surface of any plate has been 
derived by Brax and Pignol (2011) for any Ratra-Peebles 
index n. In S.I. units it reads: 



FIG. 11 Solution of the chameleon equation for n = 1 (black 
curves) and n = 2 (red curves). Solid lines correspond to the 
two-plates solution Eq. (72) with 2 R = 1 cm. Dashed lines 
correspond to the one-plate solution Eq. (71). 


ered by Ivanov et al. (2013). They proposed the approx¬ 
imate analytical solution 


<p(X) = M a 


2 + n M a \ 2/(2+n) 

V2 he J 


(71) 


<p(x) = M a 


/ 2 -(- n RM\ 
V 2y/2 he 



2/(2 +n) 


(72) 


In the particular case n — 2 we recover Eq. (56) in 
the limit ps = 0. This solution is independent of the 
coupling strength 0, because of the saturation property 
we just discussed. 

The problem of the “capacitor”, i.e. two parallel plates 
separated by a thickness 2 R of vacuum, has been consid- 


which is exact for n = 2. In this equation x = 0 refers to 
the middle position of the capacitor. 

The one-plate and two-plate chameleon profiles are dis¬ 
played in Fig. 11. What is the physical reality of those 
chameleon profiles? Consider a test particle of mass m 
in the vicinity of a plate, or inside a two-plate capaci- 

















19 


tor, the chameleon profile generates a potential energy 
(3m/m p \(p(x). We will discuss later suitable neutron ex¬ 
periments to probe these chameleon fields. 

It is important to note that a macroscopic test mass 
will not be subjected to the same potential energy. For 
strong couplings j3, the core of the test mass will be 
shielded from the external chameleon field. It is some¬ 
times useful to think of experiments searching for a new 
force as made of two components: a source (sometimes 
called an attractor) and a probe (sometimes called a de¬ 
tector). In the case of strongly coupled chameleons, a 
macroscopic source such as a plate produces a saturated 
chameleon field, independent of /3. Also, the response 
of a macroscopic probe will show some saturation effect, 
whereas the response of a particle such as a neutron will 
be proportional to the coupling strength /?. 

G. Searching for the chameleon in the lab 

Probing new forces with neutrons. Since the 
existence of new interactions are generic predictions of 
theories beyond the standard model, they are actively 
searched for in a great variety of experiments probing 
scales from subatomic to astronomical distances (see An- 
toniadis et al. (2011) for a recent review). 

It is however quite recently that slow neutrons were 
recognized as interesting probes of new forces. In par¬ 
ticular, neutron scattering experiments were found to 
be sensitive to new Yukawa forces of nanometric range 
(Nesvizhevsky et al., 2008; Pokotilovski, 2006). At the 
micrometer scale, the measurement of quantum states 
of neutrons bouncing over a mirror provides some con¬ 
straints, however they are not as stringent as those de¬ 
rived from fifth-force searches using macroscopic bodies. 

In addition, neutrons can be sensitive to spin- 
dependent forces in the sub-millimeter range (induced 
by new bosons of mass n > 10 -3 eV). Limits on 
CP violating spin-dependent forces induced by the ex¬ 
change of light Axionlike particles were first set by bounc¬ 
ing neutrons (Baefiler et al ., 2009, 2007; Jenke et al, 
2014). These limits were quickly superseded by exper¬ 
iments measuring the spin-precession of ultracold neu¬ 
trons (Afach et al., 2015; Serebrov et al., 2010). Now the 
neutron limits compete with experiments using hyper- 
polarized 3 He (Petukhov et al., 2010). At shorter range, 
a neutron diffraction experiment provides an interesting 
constraint on Axionlike particles (Voronin et al., 2009). 
There are also neutron beam experiments probing new 
spin-dependent interactions mediated by spin 1 bosons 
(Piegsa and Pignol, 2012; Yan and Snow, 2013). 

Neutrons probing strongly coupled chameleons. 

Due to the peculiarity of the force induced by the 
chameleon field, the constraints on generic short range 
Yukawa force do not generally apply. In this case, it is 


quite attractive to consider the neutron as a probe of 
strongly coupled chameleons because of the absence of 
the saturation effect. 

The case of the neutron bouncer is particularly inter¬ 
esting. Ultracold neutrons can bounce above an hori¬ 
zontal plate because they are specularly reflected by the 
Fermi potential when they fall on the plate (which we call 
a mirror). It is possible to prepare neutrons almost at 
rest (concerning the vertical motion), so that the bounc¬ 
ing height is only a fraction of a millimeter. In this case 
the energy of the vertical motion is quantized, as demon¬ 
strated experimentally for the first time at the Institut 
Laue Langevin by Nesvizhevsky et al. (2002). We (Brax 
and Pignol, 2011) argued that the chameleon field Eq. 
(71) would be detectable in this system and deduced the 
upper limit /3 < 10 11 . Next, Jenke et al. (2011) real¬ 
ized the first spectroscopy of the quantum states by in¬ 
ducing resonant transitions. The chameleon field would 
shift the frequency of the resonances and the upper limit 
/3 < 6 x 10 s was deduced (Jenke et al., 2014). In part IV 
we will dwell on this topic and give the prospects of the 
GRANIT experiment. 

Besides quantum states of bouncing neutrons, there 
is a second route: neutron interferometry. Pokotilovski 
(2013) proposed to build a Llyod’s mirror interferometer 
for very cold neutrons to probe the chameleon. Such a 
method could in principle be sensitive to chameleon cou¬ 
plings down to /3 ~ 10 7 but the cold neutron interferome¬ 
ter technique has yet to be developed. As an alternative, 
we (Brax et al, 2013) proposed to use triple Laue-case 
interferometers with slow neutrons that have been oper¬ 
ated routinely for decades in several neutron facilities. A 
first experiment has been performed in summer 2013 at 
the Institut Laue Langevin (Lemmel et al., 2015). We 
will describe it in part III. 

The constraints in the chameleon parameter space 
(Ratra-Peebles index n and matter coupling /3) from 
these neutron experiments are shown in Fig. 32 together 
with constraints from non-neutron experiments. 

Other laboratory searches for chameleons. 

Aside from neutron experiments, there are of course a 
few other means to probe the chameleon. 

• Using torsion pendulum with exquisite sensitivity 
(such as Kapner et al. (2007)) one searches for 
an anomalous torque between a rotating attractor 
disk and a torsion pendulum separated by a frac¬ 
tion of a millimeter. The detector and the attrac¬ 
tor are both disks with a diameter of a few cen¬ 
timeters, they are then subjected to the chameleon 
screening mechanisms in the case /3 1. Upadhye 

(2012) undertook the very complicated task of de¬ 
riving the limits on the chameleon couplings from 
this experiment. For n = 1 the excluded range is 
10~ 2 < /3 < 10. Couplings larger than 10 are al¬ 
lowed! 
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• Measurements of the Casimir force could also be 
sensitive to chameleons (Brax et al. , 2007). Exist¬ 
ing experiments are sensitive only in the region of 
the parameter space where M\ >0.1 eV. Future 
experiments could be sensitive to the interesting 
case M\ = 2.4 x 10“ 3 eV for a wide range of cou¬ 
pling /3. 

• Atomic precision tests in the hydrogen atom can re¬ 
veal the presence of a scalar field coupled to protons 
and electrons. Brax and Burrage (2011) derived the 
limit (3 < 10 14 from existing data. 

• An atom-interferometry experiment using ultracold 
cesium atoms (Hamilton et al. , 2015) has been re¬ 
ported (while this manuscript was nearly finished). 
They obtained the limit /3 < 2 x 10 4 for n = 1. 

• In addition to the coupling to matter (the term 
pm/mpiipijjip) , the chameleon field can be coupled 
to photons with the Lagrange density term 

l ^Lg>F flv F^. (73) 

4 TOpi 

This case leads to a unique experimental signa¬ 
ture. Due to the 77 tp vertex, chameleons can be 
produced by shining an intense laser in a region 
permeated by an intense magnetic field (take a 
spare dipole magnet from a giant particle acceler¬ 
ator). In this context, a chameleon is an elemen¬ 
tary quantum excitation of the field ip. The pro¬ 
duction region is a vacuum chamber, the produced 
chameleon has a very low mass. The walls of the 
vacuum chamber act as a repulsive potential for 
the chameleon, because the mass of the chameleon 
is large inside the walls. Therefore the chameleons 
can be stored in the chamber (provided that the 
coupling to matter is strong enough, in practice it 
requires /3 > 10 4 ). After switching off the laser, the 
trapped chameleons can be converted back to pho¬ 
tons through the reverse of the process that formed 
them, one would then measure an afterglow. A ded¬ 
icated experiment called CHASE (Chameleon Af¬ 
terglow Search) was performed at Fermilab Steffen 
et al. (2010). No afterglow was observed, the ex¬ 
periment sets the limit /? 7 < 10 11 . 

Modified gravity. We conclude this part by men¬ 
tioning that the chameleon theory is one amongst several 
known (theoretically) screening mechanisms (see Khoury 
(2013) for a pedagogical overview). In the general frame¬ 
work of modified gravity, the Einstein-Hilbert action of 
general relativity is modified by adding a scalar field p. 
There are many possible choices for the potential V(p), 
for the coupling to the standard model A(</?), for non¬ 
standard kinetic terms, etc. To understand the differ¬ 
ent possibilities we follow Brax (2013) and write the La- 



FIG. 12 Number of publications per year containing in 
the title the keywords “dark energy” (black), “cosmological 
constant” (green), “quintessence” (red), “chameleon” (blue), 
“galileon” (pink). The counting was done using inspire- 
hep. net. 

grange density linearised around a background configu¬ 
ration 

c = Z< ^-d pd^p - ^-p 2 (p 0 )p 2 + ^^-pp. (74) 

In general the background configuration p 0 depends on 
the environment density p, a property which could bring 
about screening. We can see that there are three general 
classes of screening mechanisms. 

• The mass p(po) becomes large in a dense environ¬ 
ment. This is the chameleon mechanism. 

• Z(p 0 ) becomes large in a dense environment. A 
special case of this is the Galileon field. 

• The coupling P{po) becomes small in a dense en¬ 
vironment. The Damour-Polyakov mechanism be¬ 
longs to this class. The symmetron is another ex¬ 
ample. 


The comprehensive study of these theories, their cos¬ 
mological implications, and the design of laboratory ex¬ 
periments is a vivid ongoing field of research (see Fig 
12). The hunt is on for the non-gravitational interaction 
of dark energy. 
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FIG. 13 General principle of neutron interferometry. 

III. NEUTRON INTERFEROMETRY CONSTRAINS THE 
CHAMELEON 

This part describes a search for the chameleon using 
neutron interferometry. The experiment has been re¬ 
ported also in Lemmel et al. (2015). 

A. Neutron interferometry 

The first neutron interferometer, analogous to a Mach- 
Zehnder interferometer for light, has been realized by 
Rauch et al. (1974). The principle of neutron interferom¬ 
etry is depicted in Fig. 13. A monochromatic neutron 
beam with a wavelength of A = 0.27 nm (or wavenumber 
of k = 2 tt/X = 23 nm -1 ) is split by Bragg diffraction into 
two coherent beams using a single-crystal silicon plate. 
Then part of these two beams are recombined using two 
additional similar parallel plates. The neutron detectors 
measure the flux resulting from the interference of neu¬ 
trons going through path I and path II. A phase flag 
(usually an aluminium plate with variable angle) is in¬ 
troduced in path II to record the interferogram, i.e. the 
neutron flux in detectors H and O as a function of the 
phase flag position: 

Fo = Ao + Bo cos(£fl ag ); F H = A H - B H cos(£fl ag ), 

(75) 

where £fl ag is a function of the phase flag position. Then, 
a sample is introduced in path I and a new interferogram 
is recorded: 

Fo = Ao+Bo cos(£fl ag -t-£); F H = A H —B H cos(£fl ag +£). 

, (76) 

The phase £ due to the sample is expressed as an inte¬ 
gral of the potential V ( x ) along the neutron beam in the 
sample (Rauch and Werner, 2000): 

i= -wS nx)ix - <77) 

Recording interferograms with and without sample one 
can extract the phase shift due to the sample. 


In order to probe the chameleon, we (Brax et al., 2013) 
proposed to use a rectangular vacuum cell as a sample. 
When evacuated, a chameleon field <p with a bubble-like 
profile will appear in the cell. Since the neutron has an 
energy V = firri/m v \ ip in the chameleon held, the phase 
associated with the chameleon bubble reads: 

v{x)ix (78) 

The cell needs to be evacuated in order to build up 
a chameleon held “bubble” in the cell, because even air 
at atmospheric pressure is too dense. Indeed, from the 
calculations presented in the previous part of this docu¬ 
ment, we can estimate that the chameleon held is in the 
saturated regime in a few cm of air for strong coupling 
(/3 1). The bubble solution (72) cannot build up in 

air, instead, the solution would be nearly uniform in the 
cell. Note that the neutron interferometers are generally 
operated in air. The search for the chameleon held calls 
for a dedicated experiment with a vacuum cell. 

Interestingly, it is possible in principle to switch off 
and on the chameleon bubble by adding or removing a 
small amount of gas in the cell. We can proht from this 
effect in an experiment, because this is a rather unique 
signature which is easy to exploit in a practical setup. 
If /3 = 10 s , a pressure of 10 -2 mbar of helium (about a 
million times less dense than air!) is enough to suppress 
the bubble, as we will calculate later. 

B. The experiment 

A dedicated experiment to search for the chameleon 
held was performed at the S18 instrument (see Kroupa 
et al. (2000) for a description) at ILL in summer 2013. 
We used the largest single-crystal interferometer avail¬ 
able (see picture Fig. 14) with a beam separation of 
50 mm. 

A vacuum chamber specially designed for the 
chameleon search was built by the Atominstitut group 
in Vienna. The device, depicted in Fig. 15 consists of 
an aluminium vacuum chamber with inner dimensions 
40 x 40 x 94 mm 3 and two air chambers sitting alongside 
the vacuum chamber. The whole box can be moved side¬ 
wards for swapping the vacuum cell between beam path 
I and beam path II and to probe different beam trajec¬ 
tories within the vacuum cell. The vacuum chamber is 
connected to a gas system through a hole of diameter 
5 mm on the top of the cell. The gas system comprises 
a pressure gauge, a motorized leak valve connected to 
a helium bottle and pumps. A turbomolecular pump is 
running continuously while a controlled amount of he¬ 
lium is let in through the leak valve. This way, the pres¬ 
sure of helium in the cell can be controlled in the range 
10~ 4 mbar < P < 10 2 mbar. To get higher pressure the 
turbomolecular pump needs to be disconnected. 
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FIG. 14 Picture of the large interferometer made out of a 
single-crystal silicon used in the S18 experiment. 


H detector 



FIG. 15 Top view of the setup. The chamber box (blue) can 
be moved by a remote controlled robot, allowing the beams 
to pass at different positions labeled by ’a’ to ’n’. 

During data-taking, we recorded many interferograms 
while varying the parameters of the cell (helium pressure 
and beam position). It takes about half an hour to record 
one interferogram. A major difficulty concerns the phase 
drifts due to environmental factors such as temperature 
gradients, air flow and vibrations. The use of such a 
big single-crystal interferometer is very delicate since it 
is more sensitive to these nuisances compared to smaller 
interferometers. To compensate for environmental phase 



FIG. 16 Interferograms recorded in position ’d’. The plots 
show the intensity in the detector O as a function of the path 
length difference created by the rotating phase flag. 


drifts, interferograms were recorded in parallel. It means 
that the phase flag was rotated to the first angular po¬ 
sition and neutrons were counted for a certain amount 
of time for each parameter setting. Then the phase flag 
is rotated to the next position and neutrons are counted 
again for all parameter settings etc. We recorded data 
with two different modes: the “profile mode” and the 
“pressure mode”. 

Pressure mode. In the pressure mode, the beam 
II passes in the middle of the cell (position ’d’). We 
recorded four interferograms in parallel with four dif¬ 
ferent helium pressures (2.4 x 10 -4 , 7.1 x 10 -4 , 2.7 x 
10~ 3 , 1.1 x 10 -2 mbar). A subset of these interferograms 
is shown in Fig. 16. 

Profile mode. In the profile mode, a fixed value of 
the helium pressure is established in the cell. Interfero¬ 
gram with different positions of the cell are recorded in 
parallel, to look for the bubble profile. More specifically, 
we recorded interferograms at positions ’a’, ’d’ and ’g’, 
that we will refer later as center position (y = 0, z = 
—4 mm) and side position (y = ±15 mm ,z = —4 mm). 
This mode allows to probe a wider domain of pressure be¬ 
cause the pressure needs not be changed rapidly. How¬ 
ever, there is a possible systematic effect in this mode 
associated with the variation of the wall thickness of the 
cell. The data in the profile mode require a position de¬ 
pendent phase correction based on a precise mapping of 
the wall thickness. 
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C. Expected signal 

To analyse the data, we need to calculate the field p 
inside the cell for different helium pressure. The field 
satisfies the chameleon equation A ip = ((/?). Since no 

exact analytical solution exists for this problem, a nu¬ 
merical scheme was employed. To simplify the problem, 
we considered the chameleon equation in two dimensions 
y, z transverse to the beam. Indeed, the longitudinal size 
of the cell (/ = 94 mm) is longer than the transverse size 
(a = 40 mm), we can assume that the field is approxi¬ 
mately uniform in the x direction. We solve for the field 
in units of M A , i.e. we introduce the dimensionless field 
F {y,z) = v(y,z)/M A - 

The numerical method adopted to calculate the solu¬ 
tion is described in Appendix A. Figure 17 shows selected 
results of such calculations. 

Expected chameleon phase-shift. From the cal¬ 
culations of the transverse profile F, we can extract the 
expected phase signal 


7^/3— f M A F dx 
kn z m p i J 

(79) 

m „ m , „ , 

yry/3 - M a F l e g. 

kn z mpi 

(80) 


We have introduced the effective length l e g to take into 
account the edge effects in the direction longitudinal to 
the beam. To estimate the effective length of the bubble 
in the longitudinal direction, we performed a 3D lattice 
calculation of the bubble. The 3D calculation is much 
more demanding in terms of CPU time as compared to 
the 2D calculation, in particular it does not permit to 
use a fine grid neither to repeat the calculation for mul¬ 
tiple change of parameters. Calculating with a grid size 
of 0.5 mm the empty bubble for n = 1 we found an ef¬ 
fective length of l e g = 84 mm (compare to the physical 
length of l = 94 mm). At the end, we found the following 
numerical relation between the expected phase shift (in 
degree) of the bubble £ exp and the 2D profile F: 

£ = -4.7 x 10" 9 deg x /3 x F. (81) 

The 2D profile F is a function of the transverse posi¬ 
tion y, z of the beam with respect to the cell, the helium 
pressure p in the cell and the coupling /3. 

Nuclear phase-shift. There is in principle an ex¬ 
pected phase shift proportional to the pressure due to the 
Fermi potential of the helium. Combining (2) and (77) 
we find that this phase amounts to 

P 

£nuci — A ^ 5 npie — 0.1 deg X - — , (82) 

1 mbar 



y [mm] 



y [mm] 



y [mm] 


FIG. 17 Calculations of the chameleon bubble profile trans¬ 
verse to the neutron beam for different helium pressure in the 
cell. We plot F(y, z) = <p(y, z)/Ma, for n = 1 and /3 = 5 x 10'. 
Top: pressure = 10 -J mbar, middle: pressure = 10 -4 mbar, 
bottom: in vacuum. 


where A = 0.27 nm is the neutron wavelength, l = 94 mm 
is the longitudinal length of the chamber, b = 3.26 fm 
is the bound neutron scattering length of 4 He and ?rjje 
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FIG. 18 Results of the experiment at S18. Top: phase differ¬ 
ence between center and side position in the cell as a function 
of the helium pressure (profile mode). Bottom: phase differ¬ 
ence between low pressure and 0.011 mbar (pressure mode). 
In both cases the orange line represents the theoretical expec¬ 
tation for n = 1 and j3 = 5 x 10'. 

is the number density of atoms. For pressures below 
10~ 2 mbar, as in the actual experiment, this phase shift 
is too small to be detected. We will ignore this effect. 

D. The results 

The final results of the S18 experiment are shown in 
Fig. 18. In the profile mode, there is no significant phase 
difference between the center and side position of the 
beam, even when the cell is evacuated to 10~ 4 mbar of 
helium. The data does not show the bubble profile as¬ 
sociated with a chameleon field. In the pressure mode, 
where the beam always probes the center of the cell, there 
is no significant phase difference when the helium pres¬ 
sure changes. The data is in agreement with the “no 
chameleon” hypothesis. 

To derive a limit on the coupling /?, we calculate the 
X 2 : 

x 2 m = E {m ' a 2 c,)2 < 83 > 

i i 

where the sum goes over all 8 data points Q ± Oi shown in 



FIG. 19 Calculated y 2 (/3) (according to Eq. (83)) for n = 
1,2, 3,4. The horizontal orange dashed line corresponds to 
the limit at a confidence level of 95 %. 


TABLE IV Limits on the chameleon coupling j3 (at 95 % 
C. L.) derived from the S18 experiment at ILL. 


n 

Aim 


i 

1.6 

X 

10 7 

2 

4.6 

X 

10 7 

3 

1.2 

X 

10 8 

4 

CO 

X 

10 s 


Fig. 18 thus combining the profile and pressure modes. 
Figure 19 shows x 2 (/3) calculated for n = 1,2,3,4. The 
minimum of X 2 {P) is always obtained for /3 = 0. To 
derive the exclusion limits (at a confidence level of 95 %) 
we use the standard criterion 

X 2 (Aim) = xLn + 1-96 2 , (84) 

the limits are reported in Table IV and also in Fig. 32. 
For n = 1, the limit is already a factor of 30 better than 
the current limit obtained with Gravity Resonance Spec¬ 
troscopy (Jenke et al . , 2014). 


A new experiment with the interferometer is scheduled 
for summer 2015. By improving the contrast of the in¬ 
terferometer and increasing the statistics, we can gain at 
least a factor of 10 in sensitivity. 
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IV. BOUNCING NEUTRONS WITH GRANIT 
A. Quantum states of bouncing neutrons 

If the kinetic energy of an incident neutron is less than 
the Fermi potential of a material wall (or floor, or ceil¬ 
ing) then the neutron will be reflected off the surface 
for any angle of incidence. This is the definition of an 
ultracold neutron. For example, glass has a Fermi po¬ 
tential of 90 neV. It reflects neutrons up to an energy of 
E = 90 neV, or a velocity of 4 m/s. The reflection is 
almost perfectly specular, i.e. mirror-like. The probabil¬ 
ity of non-specular (diffuse) reflection off a well-polished 
glass mirror can be less than 10 -3 (Nesvizhevsky et al . , 
2007). 

The Fermi potential of glass corresponds to a height in 
the gravity field of h = 88 cm through the relation E = 
mgh , where m is the neutron mass and g = 9.806 m/s 2 
is the acceleration in the lab. It means that a neutron 
dropped from a height of less than a meter will bounce 
above a glass floor. The bouncing period T = ^8h/g is 
of the order of a second, like any object bouncing with a 
height of about a meter. 

Now, if we consider a shorter bouncing height, or a 
higher bouncing frequency, the simple classical descrip¬ 
tion of the bouncing ball might fail. We will see that 
when the bouncing frequency enters the audio range 
(above about 100 Hz), corresponding to a bouncing 
height of about 10 pan, a quantum description is needed. 

The quantum bouncer. The vertical motion of 
a neutron bouncing above a mirror nicely realizes the 
academic problem of a particle confined in a potential 
well. The well consists in the gravitational potential 
V ( z ) = mgz pulling the particle down and a perfect mir¬ 
ror H(0) = Too 4 * * pushing the particle up. According 
to the quantum description, the vertical motion has a 
discrete energy spectrum. The bound states energies E k 
and the corresponding wavefunctions ip k {z) satisfy the 
stationary Schrodinger equation 

H 2 d 2 

- + rn.gzipk = E k ip k , if> k ( 0) = 0. (85) 

The problem of finding the bound states energies and 
wavefunctions is solvable and we are going to solve it. By 
dimensional analysis we define the characteristic height 
of the quantum bouncer: 

/ n 2 \ 1/a 

* = (*35) = 5 ’ 87 (86) 


4 The mirror potential is the Fermi potential of glass (ss 100 neV), 

it can be considered infinite because it is much larger than the 

energies of the quantum levels (« 1 peV). 


PROBABILITY FUNCTIONS IN COORDINATE SPACE 
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FIG. 20 Representation of the probability density functions 
\ip k (z)\ 2 ,k = 1,2,3, 4 for the first four energy states of the 
neutron bouncer. The four functions are translated in the 
vertical axis according to the energy of the state. The grey 
zone is classically forbidden. 

Defining the dimensionless height Z = z / zq , the station¬ 
ary Schrodinger equation takes the form of the Airy equa¬ 
tion 

+ (<* - Z)A = o, ( 87 ) 

where e k = E k /mgzQ. The general solution of this sec¬ 
ond order linear differential equation is a linear combi¬ 
nation of the two Airy functions Ai and Bi: i/j{Z) = 
C Aa(Z — e) T D Bi (Z — e). The function Bi diverges 
at Z = Too therefore the physically acceptable solutions 
are those with D = 0. Next, the condition ijj k ( 0) = 0 
at the boundary enforces the quantification of the energy 
levels Ai(— e k ) = 0. Therefore the energy levels of the 
stationary states are given by 

E k = mgz 0 e k (88) 

where 

e k = {2.338,4.088,5.521,6.787, • • ■ } (89) 

is the sequence of the negative zeros of the Airy function. 
Numerically the energy levels are in the range of 1 peV = 
10 ~ 12 eV (E 1 = 1.41 peV, E 2 = 2.46 peV •••), five 
orders of magnitude below the kinetic energy of ultracold 
neutrons. 

To each quantum state |fc) corresponds the wavefunc- 
tion 


i>k{z) = C k Ai (z/z Q - e fc ), 


(90) 









26 


where Ck is a normalization coefficient. The probabil¬ 
ity density functions | {z\ k) \ 2 = \ipk(z )\ 2 associated with 
the wavefunctions are plotted in Fig. 20. As in any ID 
potential well, the number of nodes (the places where the 
wavefunction vanishes) of the k th quantum state is k. 

Let us compare the neutron quantum bouncer to the 
hydrogen atom, a more familiar case of a bound system. 
The typical energy of the neutron bouncer is 13 orders 
of magnitude smaller than 13.6 eV, the binding energy 
of the electron in the ground state. The typical size of 
the wavefunction is 5 orders of magnitude larger than 
0.05 nm, the size of the hydrogen atom. These very un¬ 
usual scales arise due to the extraordinary weakness of 
gravity as compared to the Coulomb force. Another im¬ 
portant difference concerns the stability of the excited 
levels. An excited state of the hydrogen atom decays 
by spontaneous emission of a photon with a time scale 
of the order of 1 ns. The neutron bouncer could also 
in principle decay by spontaneous emission of gravitons. 
However, Pignol et al. (2007) estimated the timescale of 
this effect to 10' 7 s, i.e. 60 orders of magnitude longer 
than the age of the Universe. Again, this is due to the 
extreme weakness of gravity. 

Observing the quantum states. It is usually 
believed that the weakness of gravity is not the experi¬ 
mentalist’s friend. In the case of the neutron bouncer, 
however, it is a great asset. It makes the wavefunctions 
very large, almost visible to the naked eye. This property 
has been used in a series of experiments started in 1999 
at the PF2 ultracold neutron source of the ILL reactor 
(Nesvizhevsky et al ., 2003, 2005, 2002). 

The setup which served to observe for the first time 
the quantum states of the neutron bouncer is sketched 
in Fig. 21. It is a flow-through experiment: ultracold 
neutrons with horizontal velocity are coming from the 
left, they bounce on top of a horizontal, well-polished 
glass mirror. They are then counted in a gaseous detec¬ 
tor filled with 3 He. A flat neutron absorber is installed 
on top of the horizontal mirror, parallel to it. It forms 
a slit, with a height h that can be precisely regulated 
with a micrometer accuracy. The experiment consists in 
recording the neutron flux passing through as a function 
of the height of the slit. The result obtained in this ex¬ 
periment is shown in Fig. 21. When the slit height is 
h < 15 pm, no neutron passes through the slit, contrary 
to the classical expectation. In the quantum picture, the 
neutron wavefunction is a linear combination of the sta¬ 
tionary quantum states. A given state “decays” in the 
slit if its wavefunction overlaps significantly with the ab¬ 
sorber. When the absorber is lower than 15 pm, every 
state, including the ground state 11), overlaps with the 
absorber and the slit becomes completely opaque to neu¬ 
trons. This is precisely what was observed. 

A quantitative analysis of the transmission curve al¬ 
lows to extract the spatial parameters of the quantum 




FIG. 21 Top: scheme of the method used to discover the 
quantum states of the neutron bouncer. Bottom: Figure 
taken from Nesvizhevsky et al. (2002), reporting the flux in 
the detector as a function of the height of the slit. The solid 
curve is the classical fit to these data. 

states. The difficulty of the analysis lies in the choice 
of a model for the absorber. Nesvizhevsky et al. (2005) 
parametrized the transmission curve by introducing the 
critical height of the first quantum state hi = ciZq. It 
corresponds to the turning point of a classical bouncer 
with vertical energy E\. They report /q xp = (12.2 ± 
I’^syst ±0.7 stat ) pm, the systematic uncertainty takes into 
account the error in the slit size calibration and the finite 
accuracy of the absorber model. In terms of the param¬ 
eter zo, the extracted value Zq XP = (5.2 ± 0.8) pm agrees 
with the theoretical value Zq 11 = 5.87 pm. 

There is a more direct way to access the spatial struc¬ 
ture of the quantum states, by taking a photograph. In¬ 
stead of just counting the neutrons at the exit of the 
mirror, a position sensitive detector with a micrometric 
resolution could be used. As proposed by Nesvizhevsky 
et al. (2000), one should aim at observing the node of the 
second quantum state (at a height of 10 pm) which pro¬ 
vides a clear signature, instead of looking at the ground 
state which does not display a specific pattern. Selecting 
the first few levels could be done with a slit: a horizon¬ 
tal mirror and an absorber on top. At the exit of the 
slit, we place a second horizontal mirror (no absorber on 
top). The second mirror is placed 13.5 pm below the 
first mirror, in order to increase the vertical energy of 
the neutron bouncing on the second mirror. This way, 
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DIFFERENTIAL MODE 



FIG. 22 Histogram of the neutron vertical position measured 
with the nuclear track position sensitive detector. A dip cor¬ 
responding to the first node of the excited wavefuctions is 
clearly visible. The figure is taken from Pignol (2009). 


the ground state population will be suppressed and the 
population of the second quantum state enhanced. 

A first version of this experiment has been performed 
at the PF2 ultracold neutron source in 2005. A plastic 
nuclear track detector coated with a thin 235 U conver¬ 
sion layer was used. Absorption of an ultracold neutron 
in the uranium induces a fission. One of the two fission 
products is emitted in the plastic and the ionization pro¬ 
duces a latent track. After removing the conversion layer 
and etching the plastic with NaOH, the tracks become 
visible with an optical microscope. The result of this 
measurement is shown in Fig. 22. A fit of the data al¬ 
lows to extract the populations of the quantum states, 
the spatial resolution of the detector, the background 
level, the height offset and the Zq parameter (Pignol, 
2009). The result is z 0 = (6.0 ± 0.2) pm, it is more 
accurate than the transmission-through-the-slit method. 
This first measurement was not ideal because the pop¬ 
ulations of quantum states up to 110) was quite high. 
Improvements are definitely possible, by refining the pu¬ 
rity of the selection of the state |2). We aim at taking 
better pictures with the GRANIT instrument. 

Gravity Resonance Spectroscopy. Imagine we 
prepare an initial quantum state |fc). It is possible to 
induce a transition k —> l to the state |Z) with a periodic 
excitation of frequency fki = ( E — Ei)/2irh. We define 5 


fo = 


mgzp 

2ttH 


145 Hz, 


the transition frequencies write 


(91) 


fki = /o(e/c — £;)• (92) 


5 Notice that a classical ball dropped from the height zo bounces 
with a frequency of 7r X fg. 


For the first few quantum states, the transition frequen¬ 
cies lie in the audio-range, for example we expect a res¬ 
onance associated to the 1 —>■ 2 transition at the fre¬ 
quency /21 = 254 Hz. Conveniently, the audio frequency 
range is easily accessible by electrical as well as mechan¬ 
ical oscillators. This allows direct experimental access to 
the energy spectrum of the neutron quantum bouncer, a 
technique called Gravity Resonance Spectroscopy (GRS). 

Let us recall the quantum mechanics of a Rabi reso¬ 
nance, at the basis of the GRS. Assume we apply the 
external harmonic perturbation V(t) = V cos(cut), with 
the excitation frequency uj/ 27t close to a transition fre¬ 
quency fu ■ It means that the detuning Sui = ui — 2-nfki 
should be small compared to the angular frequencies of 
the transitions in the spectrum. Assume also that the 
bouncer is prepared initially in the sate |^>(0)) = |&)- 
In this case, the two level approximation holds and the 
system can be described by \4>{t)) = a{t) \k) + b(t) |Z). 
Solving the time-dependent Schrodinger equation 

ih jt = (lb +m92 ) m *)) + ^ elult ( 93 ) 


we deduce the quantum amplitudes aft) and b(t). Finally 
the transition probability P*, _>;(£) = \b(t)\ 2 is given by the 
Rabi formula: 


Pfc-w(t) 


sin 2 (\/Suj 2 + fl 2 t/2) 

1 + Su} 2 /n 2 


(94) 


where 

n=^-(l\V\k) (95) 

is the Rabi angular frequency that characterizes the 
strength of the excitation. The duration t available for 
the excitation is an important parameter. A longer du¬ 
ration is desirable because it corresponds to a narrower 
resonance line. In flow through mode, neutrons pass a 
transition region of typically 16 cm at a speed of 4 m/s, 
corresponding to a duration of t 0 = 40 ms (these pa¬ 
rameters are those of the system built for the first stage 
of GRANIT). With the duration of the excitation fixed, 
one should set the excitation strength such that Cttg = ir. 
Therefore, the transition probability is maximum at res¬ 
onance: Pfc_>z(fo) = 1 when Su> = 0. The measurement 
consists in recording the presence of the state |/) (or the 
disappearance of the state | k)) for different values of the 
excitation frequency u}/2n. It forms a resonance curve 
with a peak at lu/2it = fki . The width (FWHM) of the 
peak is A/ = fl/ir = 1/to- 

There are two practical possibilities to do the exci¬ 
tation (i) using a vibrating mirror (ii) using a magnetic 
force. GRS with vibrating mirror has been realized at the 
PF2 UCN beamline by Jenke et al. (2014, 2011). The 
technique with a magnetic excitation is in preparation 
with the GRANIT instrument. We will come back to 
this subject in detail later. 
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FIG. 23 Inertial mass versus gravitational mass of the neu¬ 
tron. The diagonal line corresponds to the weak equivalence 
principle rrii = m g . The pink band represents the la region 
deduced from the analysis of the photograph of the quantum 
states. The blue dashed line is an iso-I?o line. 

Testing the equivalence principle. We learn at 
school that it is impossible to tell the mass of a body by 
observing how it falls down. Newton’s law of motion for 
a falling body is rriia = m g g , where a is the acceleration 
of the body, g is the gravity strength, m,; is the inertial 
mass, m g is the gravitational mass. If we assume that 
both masses are equal then we find that the acceleration 
of all bodies is the same, equal to g. The hypothesis 
(rrii = m g ) is called the weak equivalence principle, 
the observable consequence is called the universality 
of free fall. 

Surprisingly, the weak equivalence principle does not 
imply the universality of free fall in the case of the quan¬ 
tum bouncer. This is because the mass in the Schrodinger 
equation (85) does not simplify, contrary to Newton’s 
law of motion. By measuring the properties of the quan¬ 
tum bouncer we can tell the mass. Indeed, from Eq. 
(86) we derive to = h/ yj2gz\. Using the result of 
the photograph Fig. 22 we measure the neutron mass: 
me 2 = 910 ± 50 MeV. 

It is therefore of interest to investigate the question 
of the equivalence principle with the quantum bouncer 
in more details. Distinguishing the inertial and gravita¬ 
tional masses, the Schrodinger equation becomes 

H 2 d 2 

- + rn g gz4>k = E k ipk, 4>k{ 0) = 0. (96) 

In this framework, the size of the wavefunctions is gov¬ 
erned by the parameter 

( h 2 \ 1/3 


It means that the measurements of the spatial properties 
of the quantum states are sensitive to the product rriim g . 
Figure 23 interprets the result of the first photograph in 
the TO, : , m g plane. Next, the energies of the quantum 
levels are E k = Eoe k with 

E ° = rn g gz 0 = ) ■ ( 98 ) 

It is remarkable that the measurement of the energy lev¬ 
els is sensitive to another combination of the masses, 
namely m 2 /rrii. 

Then, not only can we tell the mass of the neutron 
by observing the quantum bouncer, but we can measure 
separately the inertial mass rrii an d gravitational mass 
m g . For this, we need to measure the spatial feature of 
the quantum states to extract Zq on the one hand and 
to perform the spectroscopy of the quantum states with 
GRS to extract Eq on the other hand. 

B. Effects of the chameleon field 

Since the energy difference between the first quantum 
states is as low as about 1 peV, the neutron quantum 
bouncer is a sensitive system to search for hypothetical 
new forces acting on the neutron at the vicinity of the 
mirror. Such a new force could be induced by the pres¬ 
ence of the chameleon field. Since this force is attractive, 
it would shrink the wavefunctions of the quantum states 
and enlarge the energy spectrum. 

We have argued in Brax and Pignol (2011) that ex¬ 
tremely large coupling p > 10 11 would lead to new 
bound states at a distance of less than 2 pm which is 
ruled out by the first transmission-through-the-slit ex¬ 
periment. The limit p < 10 11 is practically independent 
of the Ratra-Peebles index n. Next, Jenke et al. (2014) 
analyzed their experiment with the vibrating mirror to 
derive limits on the chameleon coupling. They quote the 
limit /3 > 6 x 10 s , also practically independent of n. 

Let us now analyse quantitatively how the properties 
of the quantum bouncer are modified in the presence of 
the chameleon field. As we have seen in part II, if the 
chameleon is strongly coupled with matter (/3 1), a 

scalar field ip(z) will be generated above the mirror. We 
have seen that the field ip(z) is saturated: it is indepen¬ 
dent of P as long as P 1, it is also independent of the 
mass density of the mirror. When no absorber is placed 
above the mirror the scalar field is given by Eq. (71). 
Then, the neutron feels the potential 

Tfi 

V(z) = mgz + p - <p(z). ( 99 ) 

TO.pi 

The chameleon term in the potential has the form 

777 

SV(z) = P - ip(z) = pv n {z/\) a 


z 0 = 


2 mim g g 


( 97 ) 


( 100 ) 
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Shifts of the transition frequencies. At first 
order in perturbation theory, the correction of the energy 
levels is 


5E k = /3v n (z 0 /X) an O kk (a n ), (104) 


and the correction of the transition frequencies is 


Sfki 


SE k - SEi 
2ttH 

j3^E(z 0 /\) a -(O kk (a n ) - 0„(«„)). (105) 


For example, consider the shift of the resonance line 2 —> 
1 (expected at /21 = 254 Hz) due to the chameleon with 
the Ratra-Peebles index n = 2. We find numerically 


Sf 21 =/3x 3.7 x 1CT 8 Hz. (106) 


FIG. 24 Wavefunction of the second quantum state ■02 ( 2 ) ■ 
The black curve corresponds to the unperturbed case (no 
chameleon). The red curve is the wavefunction in the pres¬ 
ence of the chameleon with n = 2 and j3 = 10 9 , calculated at 
the first order of perturbation theory. 


with 


v n 


m 


M a 


m pi 



2 + n 


Oin 

x 1CT 21 eV, 
( 101 ) 


where A = Kc/M\ = 82 pm and a n = 2/(2 + n). 

The task is now to calculate the shift of the energy 
states and the shrinking of the wavefunctions induced by 
the perturbation 5V. We skip the details of the calcula¬ 
tions (see appendix B) and we discuss the results. 


Shrinking of the wavefunctions. At first order in 
perturbation theory, the correction of the wavefunctions 
is 


i>k(z) = ^ k \ z ) +/3v n (z 0 /X) a " Y 


Okl(an) 
E k — Ei 



(z). 


( 102 ) 

In Fig. 24, we plot the modification of the second state 
for the coupling (3 = 10 9 using the results in the appendix 
B (the sum goes up to l = 5). The effect of the chameleon 
is a shrinking of the wavefunction, as expected for an 
attractive force. The node of the wavefunction hi (which 
we will try to measure in GRANIT) is displaced to lower 
2 . 


More quantitatively, a numerical study in the case n = 
2 shows that the displacement of the node is a linear 
function of (3 for (3 < 10 9 : 


h° 2 = (10.3 - /3 x 7.7 x 1CT 10 ) pm (103) 


A precision of 0.1 pm, which is what we could reason¬ 
ably get in forthcoming measurements with GRANIT, 
will correspond to /? ~ 10 s . This method is less sensitive 
than the gravity resonance spectroscopy, as we will see. 


As we will see later, we aim at a precision better 
than 1 Hz for the measurement of the 2 —> 1 transi¬ 
tion frequency with the flow-through resonance setup in 
GRANIT. This measurement will therefore be sensitive 
to a chameleon coupling of j3 « 2 x 10'. 

It is possible in principle to improve the precision on 
the transition frequency by increasing the interaction 
time t. We recall that the width of the resonance is given 
by A/ = 1/f, and the precision is approximately given 
by = AZ/v^, where N is the number of counts in 
the resonance. Instead of a flow-through setup, one needs 
to trap the neutrons in quantum states. Ultimately, the 
precision is limited by the beta decay lifetime of the neu¬ 
tron T n ss 880 s. Therefore the ultimate precision is 
Sf S3 1 /(r n y/N) ss 10 -5 Hz (taking N = 10000). In con¬ 
clusion, the measurement of the spectrum of the quan¬ 
tum states could potentially be sensitive to chameleon 
couplings down to (3 = 10 2 (a long way to go, though...). 


C. The GRANIT instrument 

GRANIT is an instrument dedicated to the study of 
the neutron quantum bouncer, now in the commissioning 
phase. For a description of the instrument see Baefiler 
et al. (2011a,b); Roulier (2015); Roulier et al. (2015); and 
Schmidt-Wellenburg et al. (2009). GRANIT is located at 
the level C of the ILL reactor in Grenoble. As shown in 
Fig. 25, the instrument comprises a superthermal 4 He 
source (fed by the monochromatic cold neutron beam 
H172A), producing ultracold neutrons to be used in the 
spectrometer located in a ISO 5 class clean room. 

No attempt to describe the instrument in technical de¬ 
tails will be made here, except for a few selected items. 
We will first describe the ultracold neutron source. Next, 
the setup for the magnetic excitation of resonant transi¬ 
tions will be explained. 

The ultracold neutron source. Previous experi¬ 
ments observing the neutron quantum bouncer were all 
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FIG. 25 Top: schematics of the GRANIT instrument at the 
ILL H172A beam line. Bottom: picture of the spectrometer 
inside the clean room. 


performed at the PF2 instrument, the “first generation” 
ultracold neutron source situated at the level D of the 
ILL reactor. “Next generation” sources can potentially 
deliver many more ultracold neutrons, and GRANIT op¬ 
erates now a source based on the SUN1 prototype de¬ 
scribed in Piegsa et al. (2014) and Zimmer et al. (2011). 

The concept of the source has been proposed by Golub 
and Pendlebury (1977): a volume filled with superfluid 
4 He is irradiated with a beam of cold neutrons. Neu¬ 
trons can undergo inelastic scattering in the medium, 
losing energy while emitting a phonon. Given the dis¬ 
persion relation of the phonons in superfluid helium, it 
turns out that a neutron having a wavelength of exactly 
A* = 0.89 nm transfers its entire energy to the phonon 
and comes at rest in the superfluid. Neutrons with a 
wavelength not too far from A* come almost at rest. The 
process produces ultracold neutrons, with a volumic pro¬ 
duction rate (number of UCN produced per unit volume 
and unit time) of 


p= W 

dX 


x (3.0 ± 0.2) x 10 9 nm cm 1 , 


(107) 


where d(f>/d\ is the differential flux in the cold neutron 
beam. This theoretical estimate is valid for neutrons cre¬ 
ated in the velocity range 0 < v < 5.6 m/s (Piegsa et al . , 
2014), corresponding to the storage capability of a trap 
made with stainless steel walls 6 . 

When the reactor power is 58 MW, H172A deliv¬ 
ers 3 x 10 9 neutrons per second with a wavelength of 
0.89 nm (the differential flux has a Gaussian shape 
with a RMS a = 0.021(1) nm), in a square section of 
7x7 cm 2 (Roulier, 2015). In the continuity of the 
guide, the source itself consists in a tube with square 
section (inner dimensions 7 x 7 x 100 cm 3 ) with BeO 
walls filled with superfluid helium. The differential flux 
inside this production volume is reduced by absorption 
of cold neutrons in windows (x0.9) and by the diver¬ 
gence of the beam (x0.8), resulting in an effective flux 
of d(j)/dX\ e tt = 8 x 10 8 cm _2 s _1 nm _1 . The UCN produc¬ 
tion inside our source is then expected to occur at a rate 
of P « 2.3 UCN/cm 3 /s. In total the 4900 cm 3 should 
produce about 10,000 ultracold neutrons per second. 

Those ultracold neutrons produced in the superfluid 
helium are stored by the BeO walls of the production 
volume. To extract the neutrons, the volume is connected 
to a chimney (a stainless steel tube) going up to a UCN 
valve. When the valve is closed UCNs accumulate in the 
volume. Above the valve there is a UCN guide going up, 
then a 90 degree bent, followed by an horizontal stainless 
steel guide connected to the spectrometer inside the clean 
room. 

For efficient accumulation of UCNs inside the source, 
the storage time r should be as long as possible in order 
to maximize the UCN density p = Pt at saturation. The 
loss rate 1/r inside the source has several contributions: 

1/t = 1/t„ + l/r wa n + 1/rsHe + 1 /t up (108) 

where r n ss 880 s is the beta decay lifetime, r wa n corre¬ 
sponds to the losses at wall collision, T3 He is due to the 
absorption by residual traces of 3 He in the bath of nor¬ 
mally pure 4 He, and r up corresponds to upscattering by 
thermal excitations in superfluid helium. Inconveniently, 
this last contribution forces us to operate the source at 
temperatures below 1 K. According to the theoretical es¬ 
timation by Golub (1979), the temperature dependent 
upscattering rate is 

l/r up = Ae~ E * /T + BT 7 + CT 3 / 2 e~ A/T (109) 

where the first term corresponds to one-phonon absorp¬ 
tion with A = 500 s — , the second term corresponds to 


Our trap is made of BeO walls with a Fermi potential of 257 neV, 
but a small portion is made of stainless steel with a Fermi po¬ 
tential of 184 neV only. Since superfluid helium has a Fermi 
potential of 18 neV, the threshold energy for trapped neutrons is 
E < 184 — 18 neV, corresponding to a critical velocity of 5.6 m/s. 


A* 
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FIG. 26 Top: Theoretical estimation of the lifetime due to 
upscattering r up in superfluid helium, as a function of the 
temperature. Bottom: Flux of UCN at the measured at the 
exit of the horizontal extraction guide (a tube with inner di¬ 
ameter 3 cm) as a function of the temperature of the helium 
(with the UCN valve completely open). 


two-phonon process with B = 0.008 s _1 , the last term 
corresponds to roton-phonon process with C = 18 s -1 . 
Helium temperature T is expressed in Kelvin, E* = 12 K 
is the phonon energy corresponding to A* and A = 8.6 K. 
The time r up is plotted as a function of the tempera¬ 
ture in Fig. 26 (top). We see the dramatic effect of the 
temperature in Fig. 26 (bottom): we measured the flux 
of UCNs at the exit of the source as a function of the 
temperature. The flux is significantly reduced when the 
temperature is above 1 K. Even at low temperature we 
only count 400 UCNs per second. Some losses are due to 
absorption of neutrons in the vacuum separation window 
(a 15 pm thick titanium foil situated in the horizontal 
guide) and the entrance window of the detector (another 
15 pm thick titanium foil). Also, the storage time at low 
temperature is only about 30 s, which is not very long 
compared to the emptying time of 13 s (Roulier et al ., 
2015). There are also certainly some more losses which 
are being investigated. 
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FIG. 27 Scheme of the connection of the source to the spec¬ 
trometer. Not shown: container of the superfluid, thermal 
screens, vacuum chamber of the source, lead biological shield, 
wall of the clean room, vacuum chamber of the spectrometer, 
support of the copper antichamber and semi-diffusive slit. 


Let us describe a little bit what’s inside the spectrom¬ 
eter. As shown in Fig. 27, the horizontal guide is con¬ 
nected to a cylindrical storage volume labeled “copper 
antichamber” (in the test mentioned before, a detector 
was placed instead of the antichamber). The 30 cm long 
cylinder has a diameter of 4 cm, it is made of pure cop¬ 
per (Fermi potential 168 neV). A small aperture in the 
antichamber leads to a semi-diffusive slit made of a bot¬ 
tom reflective mirror and a top diffusing mirror. Both 
mirrors are coated with diamond like carbon. The idea 
is that neutrons with too large vertical velocity in the slit 
are diffused back in the antichamber. On the contrary, 
neutrons with small vertical velocity, those bouncing on 
the bottom mirror without touching the roof, are trans¬ 
mitted. In a first test of this system conducted in De¬ 
cember 2014, the opening of the semi-diffusive slit was 
set to 127 pm, preparing the first 20 quantum states. 
A 3 He counter with an aluminium entrance window was 
placed immediately at the exit of the slit. Unfortunately 
due to a cryogenic failure it was only possible to get to 
T = 1.3 K for the temperature of the source. In this con¬ 
figuration, we counted 0.13 neutrons per second out of 
the semi-diffusive slit. This is not yet as good as the ex¬ 
periments performed at PF2 (see Fig. 21). By fixing the 
cryogeny and using a detector with improved efficiency 
for soft UCNs (i.e. with a titanium entrance window for 
example) things will get better. Progress is under way to 
deliver a decent flux of ultracold neutrons at the exit of 
the slit. 


The magnetic excitation. Now I describe how 
we plan to measure the resonant transitions using the 
magnetic excitation as proposed in Baefiler et al. (2015); 
Kreuz et al. (2009); Pignol (2009); and Pignol et al. 
(2014). An oscillating magnetic field gradient will gener¬ 
ate the vertical force on the neutron to excite the system. 
Indeed, the motion of a neutron in a magnetic field B is 
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FIG. 28 Sketch of the flow through setup. UCNs enter from 
the semi-diffusive slit on the left ; they fall down the step to 
depopulate the ground state (1), 16 cm long transition region 
(2), 9 cm long filter (3), and detector (4). 


governed by the potential s/.t|B| 7 , where fj, = 60 neV/T 
is the magnetic moment and s = 1 for “spin up” neutrons 
and s = — 1 for “spin down” neutrons. Classically, a ver¬ 
tical force is applied on the neutron by the field gradient 
cb|B| (an homogeneous field does not exert any force, 
it acts only on the spin dynamics). Let us now assume 
a magnetic excitation of the form |B| = fdz cosuit, the 
corresponding quantum mechanical potential reads 

V(t) = s^ifdz cosut. (HO) 

This is what we need to excite resonant transitions be¬ 
tween the quantum states of the neutron bouncer. 

The sketch of the setup is shown in Fig. 28. The spec¬ 
troscopy method is conceptually the same as the famous 
1939 Rabi resonance experiment. It is performed with 
four steps (1) state preparation, (2) resonant transition, 
(3) state analysis, and (4) detection of the transmitted 
flux. 

1. UCNs are first prepared in an excited state by going 
down a step (1) of height 15 pm. The population 
of the first state |1) is suppressed as compared to 
the populations of the excited states |2), |3), • • •. 

2. Next, transitions between quantum states are ex¬ 
cited with a periodic magnetic field gradient. The 
length of the transition region L = 16 cm corre¬ 
sponds to an average passage time to = 40 ms 
8 . Two different schemes could be implemented in 
principle: the AC excitation ad the DC excitation. 
In the DC mode, the field gradient is static and spa¬ 
tially oscillating in the x direction with a period of 
d = 1 cm. In this case, only neutrons with specific 
horizontal velocities meet the resonance condition. 
The deexcitation 2 —> 1 is expected to occur at a 
frequency of / 2 1 = 254 Hz corresponding to the res¬ 
onant horizontal velocity of V 21 = d /21 = 2.54 m/s 
(in the 3 —)• 1 case we expect V 31 = 4.62 m/s). In 


' This is valid when the variation of the magnetic field is slow com¬ 
pared to the Larmor frequency. In this case, called the adiabatic 
regime, the spin follows the direction of the magnetic field and 
the neutron trajectory and the spin dynamics are decoupled. 

8 We expect the v x velocity along the beam to be distributed with 
a mean value of 4 m/s and a standard deviation of 1.5 m/s. 



FIG. 29 Picture of the wire array for the magnetic excitation. 
It is made of four modules, each one holding 32 adjacent wires 
with a section of 1 mm 2 each and a length of 30 cm in the y 
direction. Adjacent wires are separated by a gap of 0.25 mm. 


the AC mode, the field gradient is spatially uni¬ 
form and oscillating in time. One would then find 
the resonances by directly scanning the excitation 
frequency. 

3. A second horizontal mirror above the main mirror 
serves as a state filter. For a slit opening of about 
25 pm, only neutrons in the ground state are ac¬ 
cepted and higher quantum states are absorbed. 

4. Neutrons are finally detected at the exit of the fil¬ 
ter. In the AC mode one should see a resonance 
in the transmitted flux as a function of the exci¬ 
tation frequency. In the DC mode, one needs to 
measure the horizontal velocity of transmitted neu¬ 
trons. This is achieved by measuring the height of 
the neutrons after a 30 cm long free fall. 

We need to apply an oscillating gradient with an am¬ 
plitude of j3 = 0.22 T/m to induce the 2 —> 1 transition 
and /3 = 0.74 T/m for 3 —> 1 (Pignol et al., 2014). This 
is done with an array of 128 copper wires with square 
section arranged as shown in Fig. 28. The array will 
be placed above the horizontal mirror in region ( 2 ) at 
a distance of 0.8 mm from the mirror. The real system 
shown in Fig. 29. Electrical connectors are arranged 
so that the following 8 -periodic pattern current could be 
applied: Ji, I 2 , h, h, ~h, ~h, ~h, ~h, h, ■ ■ ■ and thus 
the magnetic field produced by the array will be 1 cm 
periodic in the x direction. It is possible to tune the 
currents in the 4 circuits I±, / 2 ,I 3 , I 4 to create a homo¬ 
geneous gradient at the surface of the mirror, as shown 
in Fig. 30. The configuration of the gradient depends on 
the external magnetic field: by applying a strong exter¬ 
nal field (of the order of 2 rnT) on top of the field created 
by the wire array, one generates an oscillating gradient 
in the x direction. The wire array is therefore a versatile 
device to generate the field gradient that can be used for 
the AC as well as for the DC excitation modes. 
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FIG. 30 Calculation of the magnetic field gradient <9-|B| 
produced at the surface of the main mirror with currents set 
to Ji = / 4 = 1.4 A, I 2 = I 3 = 3.5 A and without any external 
field applied. 

We will apply currents up to 5 A in the circuits, corre¬ 
sponding to Joule heating of up to 25 W (the resistance 
of each circuit is 0.3 fl). Since the wire array will be op¬ 
erated in vacuum, an active cooling system is necessary 
to dissipate the heat. We built a heat exchanger with cir¬ 
culating cold gaseous nitrogen and tested it successfully 
in vacuum. 

In the AC mode, we apply AC currents in the wires. 
Note that the gradient /3 is not a linear function of 
the currents. In particular, the gradient always satis¬ 
fies /3 > 0, because the field |B| is always weaker away 
from the wire. As a result, the excitation frequency (i.e. 
the frequency of (3(t) seen by the neutrons) is twice the 
driving frequency (i.e. the frequency of the AC current 
in the wires). The complete system of excitation has 
been simulated in Pignol et al. (2014) by solving the 
time-dependent Schrodinger equation with realistic func¬ 
tion /3(f) taking into account the non-linearities. The 
expected resonance curve is shown in fig. 31. We see 
a double resonance, one corresponding to spin up neu¬ 
trons, the other to spin down. This is a shift called 
Stern-Gerlach shift (Baefiler et al, 2015) caused by the 
unavoidable constant (time-independent) component of 
/3(f). Fortunately this shift cancels when averaging the 
two resonant frequencies. 

With a statistics of 10,000 neutrons counted during a 
scan of the driving frequency between 50 Hz and 200 Hz, 
it will be possible to extract the resonance frequency /21 
with a precision better than 1 Hz. 



FIG. 31 Calculation of the expected resonance curve in the 
AC mode of excitation. The probability of the 2 —> 1 transi¬ 
tion is plotted as a function of the driving frequency. 

the neutrons in the quantum states for long times com¬ 
pared with the flow-through arrangement. One would 
use a trap consisting of an horizontal mirror and verti¬ 
cal walls, the excitation (vibration or magnetic) could be 
applied during a longer storage period. The control of 
the environmental noises such as vibrations is an issue 
because it could provoke unwanted transition between 
quantum states, but the lifetime of the quantum states 
in the trap could reach a few tens of seconds (Codau 
et al., 2012; Pignol, 2009). The filling and emptying of 
such a trap with limited losses of neutrons remains an 
open problem. 

Another interesting route for the future is the devel¬ 
opment of position sensitive detectors with micrometric 
resolution to “see” the quantum states. We have men¬ 
tioned an attempt to take pictures of the quantum states 
using a uranium-coated plastic nuclear track detector. It 
achieves a good resolution of 2 pm, but it lacks real-time 
readout. More recently, pixel detectors based on CCD 
(charge coupling device) coated with 6 Li or 10 B have 
shown to offer both a good spatial resolution and real¬ 
time readout capability (Kawasaki et al, 2010). Pursuing 
in this way will certainly lead to clean measurements of 
the wavefunctions. 

In conclusion, progress are foreseen both in the spec¬ 
troscopy of the quantum bouncer and in the measure¬ 
ment of the spatial features of the quantum states. These 
will permit to probe new forces acting on the neutron and 
to test the equivalence principle in a quantum context. 


In the long term, increasing the precision of the tran¬ 
sition frequencies can in principle be achieved by storing 
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V. CONCLUSION 

It is fascinating that experiments at low energy with 
neutrons can address three big questions about the Uni¬ 
verse: 

1. Does the standard theory of primordial nucleosyn¬ 
thesis predict the correct amount of helium and 
deuterium? 

2. Was the asymmetry between matter and antimat¬ 
ter generated during the electroweak phase transi¬ 
tion? 

3. Is Dark Energy a dynamical field interacting with 
matter? 

The measurement of the neutron lifetime helps to an¬ 
swer question 1. Experiments will continue until the in¬ 
consistencies of the present data are resolved. 

The answer to the second question is either yes or no. 
By improving the measurement of the neutron electric 
dipole moment we could be able to reach a definitive 
answer in the next decade. 

The third question is the main focus of this work. We 
have looked in details at the strongly coupled Khoury- 
Weltman chameleon. The present constraints are com¬ 
piled in Fig. 32. We have identified two methods to probe 
the chameleon with neutrons. The first is neutron inter¬ 
ferometry. A pilot experiment performed in 2013 was 
described in this work. The second is using the quan¬ 
tum states of bouncing neutrons. Ultimately, it has the 
potential to explore a significant part of the parameter 
space. 

Searching for Dark Energy in the lab is a new and 
rapidly developing field of research. Less than ten years 
after the chameleon was proposed, a major part of the pa¬ 
rameter space is already excluded by experiments search¬ 
ing for deviations of the inverse square law of gravity at 
short distances, neutron experiments, and very recently 
atom-interferometry. Besides the chameleon, there are 
other theoretical ideas related to Dark Energy to be ex¬ 
plored in the lab. For sure there are other neutron ex¬ 
periments to be invented. 



n 


FIG. 32 Exclusion regions (95 % C.L.) in the chameleon 
parameter space (Ratra-Peebles index n and matter coupling 
1 3 ). The blue zone is excluded by gravitational resonance spec¬ 
troscopy of the quantum bouncer (Jenke et al, 2014). The red 
dashed line corresponds to the sensitivity of the flow-though 
setup in GRANIT, the red dotted line corresponds to the ulti¬ 
mate sensitivity. The orange zone is excluded by the neutron- 
interferometry experiment (Lemmel et al, 2015). The grey 
zone is excluded by the Eot-Wash experiment to search for 
a short range force (Upadhye, 2012). The hatched region is 
excluded by the atom-interferometry experiment (Hamilton 
et al, 2015). 
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Appendix A: Lattice calculation of the chameleon field 


Appendix B: Perturbation of the quantum states 


We describe here the calculation of the chameleon field 
in the rectangular cell used in the neutron-interferometry 
experiment. The dimensionless field F = <p/M\ satisfies 
the equation 

d 2 F d 2 F T ~( n Pph 3 c 5 \ 

dy 2 + dz 2 ~ ( F n + l + M 3 top J ( ^ 


We use the Rayleigh-Schodinger perturbation theory 
to calculate the effect of the chameleon field on the quan¬ 
tum states of the bouncing neutron. The perturbation 
5V is given by Eq. (100). The shift SE k of the energy 
level E k is given at first order by 

5E k = {k\ SV(z) | k) (Bl) 


with K = = (0.082 mm) -1 . The second term on 

the right hand side of the equation is proportional to the 
helium pressure (p is the mass density of helium in the 
cell). In addition, we assume that the field is saturated 
in the walls of the cell, therefore we impose F = 0 at the 
boundary of the cell. 

Let us now describe the numerical method adopted to 
calculate the solutions. We used a finite difference on a 
2D grid of step h. The Laplace operator on the grid is 
taken to be 

{AF) i j = — (F i+ i t j + Fij +1 + Fi_ij + F i: j_i — 4 Fij), 

(A2) 

then the discrete chameleon equation becomes 

4 F ii:j = Nij - h 2 Sij (A3) 


where 


Aj F.^i j -f- F^ j _|_i 


Fi 


*—i,j 


■ Fi 


i,j —1 


Sij =K 2 - 


f3p{hc) 3 \ 

F " +1 ' M A m P i J ' 


+ 


We used the iterative Jacobi relaxation method with SOR 
(Successive Over-Relaxation). The algorithm starts from 
an initial field F^ ( in our case the initial field is con¬ 
structed from the ID profile (72)) and iterates according 
to 


' l ,3 


(1 -w)fg } 


UJ 


+ 7 A, 


(fc) 
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(A4) 


The boundary condition (F t j = 0 for i,j at the bound¬ 
ary) is enforced at each iteration. If the sequence con¬ 
verges, the fixed point necessarily satisfies the discrete 
chameleon equation. The SOR parameter to usually 
serves to accelerate the convergence (in this case one sets 
u) > 1). In our case, we observed numerical instabilities 
spoiling the convergence when calculating the bubble for 
high helium pressure. The numerical instabilities were 
fixed by setting the SOR parameter to oj = 0.8, at the 
price of a slower convergence. Finally we used a step of 
h = 0.2 mm and iterate the algorithm 30000 times. A 
calculation takes about 10 min in a normal computer. 
Figure 17 shows selected results of such calculations. 


where | k) are the unperturbed stationary states. Like¬ 
wise, the modification of the wavefunction at first order 
of the perturbation theory is 


i’k(z) = + 

l^k 


(l\ SV(z) \k) 

—S-p—W (-)• 


(B2) 


where ipj.°\z) = (z\ k) are the unperturbed wavefunc- 
tions. 

It is useful to define the matrix elements of the power 
law operator as: 

/»00 

O k i(a) = (l | ( z/z 0 ) a | k) = / i>i(z)ip k (z)(z/z 0 ) a dz. 

Jo 

(B3) 

These matrix elements are calculated numerically for the 
first five quantum states and for the Ratra-Peebles in¬ 
dices n = 1, 2, 3,4. We provide the result in the form of 
four symmetric 5x5 matrices: 
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